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The 2-Channel Kondo Model II: CFT Calculation of Non-Equilibrium Conductance 
through a Nanoconstriction containing 2-Channel Kondo Impurities 

Jan von Delft^O, A. W. W. Ludwig^, Vinay Ambegaokar^ 

^Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, NY 14-853, USA 

^ University of California, Santa Barbara, CA 93106 , USA 

(September 9, 1998. Published in Ann. Phys. 273, 175-241 (1999).) 

Recent experiments by Ralph and Buhrman on zero-bias anomalies in quenched Cu nanocon- 
strictions (reviewed in a preceding paper, I) are in accord with the assumption that the interaction 
between electrons and nearly degenerate two-level systems in the constriction can be described, for 
sufficiently small voltages and temperatures {y,T < Tk), by the 2-channel Kondo (2CK) model. 
Motivated by these experiments, we introduce a generalization of the 2CK model, which we call the 
nanoconstriction 2-channel Kondo model (NTKM), that takes into account the complications arising 
from the non-equilibrium electron distribution in the nanoconstriction. We calculate the conduc- 
tance G{V, T) of the constriction in the weakly non-equilibrium regime of F, T <C Tk by combining 
concepts from Hershfield's F-operator formulation of non-equilibrium problems and Affleck and 
Ludwig's exact conformal field theory (CFT) solution of the 2CK problem (CFT technicalities will 
be discussed in a subsequent paper. III) . Finally, we extract from the conductance a universal scaling 
curve r(t') and compare it with experiment. Combining our results with those of Heftier, Kroha 
and Hershfield, we conclude that the NTKM achieves quantitative agreement with the experimental 
scaling data. 

PACS numbers: 72.15.Qm, 72.10.Fk, 63.50.-fx, 71.25.Mg 
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I. INTRODUCTION 

This JSr-, the second in a series of three papers 
(I,II,III)ErEl devoted to the 2-channel Kondo model 
(2CK). In the preceding paper (I), we gave a detailed 
review of a possible experimental realization of this 
model, jjiamely the experiments by Ralph and Buhrman 
(RB)q Q on non-magnetic zero-bias anomalies (ZBAs) in 
Cu nanoconstrictions, ancLxelated experiments by Upad- 
hyay, Louie and BuhrmanQ on Ti constrictions. The ex- 
perimental facts were summarized in the form of fifteen 
important properties of the data, nine for Cu and six 
for Ti constrictions, [see (Cu.l) to (Cu.9) and (Ti.l) to 
(Ti.6), in section IV of I). The main conclusion of paper I 
was that all experimental facts are in accord with the as- 
sumption that the ZBA is caused by the scattering of 
electrons off nearly degenerate two- level systems (TLS), 
with whom they interact acccupding to the non-magnetic 
Kondo model of Zawadowskulld, which renormalizes to 
the 2CK model at sufficiently low temperatures. (See 
Appendices y| and |Q for background on Zawadowski's 
model.) 

In the present paper (II), we focus on property (Cu.6) 
for the Cu constrictions: in the so-called weakly non- 
equilibrium or non- Fermi-liquid regime of sufhciently 
small voltages and temperatures {V <^ Vk and T <C T^, 
but arbitrary ratio v = eV/kgT) where Vk and T^ are 
experimentally determined cross-over scales) the conduc- 
tance CUV, T) was found to satisfy the following scaling 
tionQ: 



relati 



G{V,T) - G{0,T) 

'J'a 



F{v) 



(1) 



with scaling exponent a — i^- This was interpreted 
as strong evidence that the samples fall in the low- 



temperature regime of the 2CK model, because its con- 
formaljfield theory (CFT) solution by Affleck and Ludwig 
(AL)ll2lO suggests precisely such a scaling form near its 
T = fixed point, and correctly predicts that a — ^, a,s 
observed. 

If this interpretation is correct, it would imply that 
RB had directly observed non-Fermi-liquid behavior, be- 
cause in the 2CK model, the exponent a = ^ is one of 
the signatures of non-Fermi-liquid physics (for a Fermi 
liquid, a = 2). Thus RB's experiments attracted a lot of 
interest, because non-Fermi-liquid behavior, so treasured 
by theorists, has been rather difficult to demonstrate un- 
ambiguously experimentally. 

However, it is of course quite conceivable that the scal- 
ing behavior can also be accounted for by spme other 
theory. Indeed, Wingreen, Altshuler and Meirll^^"' have 
pointed out that an exponent of a = 2 ^-Iso arises within 
an alternative interpretation of the experiment, based not 
on 2CK physics but the physics of disorder (We believe, 
though, that their scenario contradicts other important 
experimental facts, see Appendix A.l of Paper I. More- 
over, the agreement between the exjperimental and the- 
oretical scaling curves found in RefJl3(a) was,the result 
of a recently discovered error in the analysislij; once the 
error is corrected, the agreement ceases.) 

It is therefore desirable to develop additional quan- 
titative criteria for comparing the experiment to var- 
ious theories. Now, in paper I it was shown that a 
sample-independent scaling function T{v) could be ex- 
tracted from the sample-dependent scaling function F{v) 
of Eq. (n[). According to the 2CK interpretation, this 
y{v) should be a universal scaling function, a fingerprint 
of the 2CK fixed point, independent of sample-specific 
details. A very stringent quantitative test of any theory 
for the RB experiment would therefore be to calculate 
r(i;), and compare it to experiment. 

The present paper is devoted to this task. T{v) is calcu- 
lated analytically within the framework of the 2CK model 
and its exact CFT solution by AL, and the results are 
compared to the RB experiment. When combined with 
recentj|Uj|merical results of Hettler, Kroha and Hershfield 
et a/Jlj'll3, agreement with the experimental scaling curve 
is obtained, thus lending further quantitative support to 
the 2CK interpretation for the Cu constrictions. 

In order to describe the scattering of electrons off 
two-level systems in a nanoconstriction geometry, we 
introduce a generalization of the 2CK model, which 
we call the nanoconstriction two-channel Kondo model 
(NTKM), that takes into account the complications aris- 
ing from the non-equilibrium electron distribution in the 
nanoconstriction. The generalization consists of labelling 
the electrons by an additional species index a = {L,R), 
which denotes their direction of incidence (toward the 
left or right for electrons injected from the right or left 
lead). 

In equilibrium {V — 0), our NTKM reduces to the 
2CK model. Therefore, for T < T^, it displays the 
same non-Fermi-liquid behavior as the latter. When 



the voltage is turned on, by continuity there must ex- 
ist a regime in which the voltage is still sufficiently small 
(namely V <ti T^) that non- Fermi-liquid behavior per- 
sists despite V ^ 0. We shall call this T,V ^T^ regime 
the non- Fermi-liquid regime, and associate it with the 
scaling regime of (Cu.6) identified in the experiment. At 
higher voltates {V ~ T^), the non- Fermi-liquid behav- 
ior is destroyed. Therefore, we shall focus exclusively 
on the case V <^ Tj^ in this paper, and accordingly the 
acronym NTKM will henceforth be understood to stand 
for "the nanconstriction 2-channel Kondo model in the 
non-Fermi-liquid regime" . 

The non-Fermi-liquid regime has to be treated by 
non-perturbative methods. The method we use com- 
bines ideas from CFT with concepts from Hcrshficldls 
y-operator formulation of non-equilibrium problemsliH. 
We show that all one needs to calculate the current us- 
ing Hershfield's formalism are certain scattering ampli- 
tudes, to be denoted by f/^,,'(e'). We assume that in 
the non-Fermi-liquid regime, the scattering amplitudes 
are essentially independent of V (since T^-dependent 
corrections are of order V/T^ <C 1 and hence neg- 
ligible (they are discussed in Appendix |). We then 
show that the V = values of the scattering ampli- 
tudes can be extracted from an equilibrium Green's func- 
tion Gjpj'{T,ix;T' ,ix') — —{T^p,i{T,ix)^p,{T,ix')) that is 

known exactly from CFT. Once the C/^^'(e') are known, 
it is straightforward to calculate the non-linear current 
I{V,T) through the constriction, and extract from it the 
scaling function T{v). 

No knowledge of CFT is required to read the present 
paper (excepting Appendix |), because the only step for 
which CFT is really needed, namely the calculation of 
Griri' , is carried out in paper III, and here we only cite the 
needed results. (Actually, the calculation of G^,,', too, 
can be done without CFT, since it has been shown very 
recently using abelian bosonization that AL'«-|&een's 
functions can be obtained without using CFTll3't2l; this 
will be discussed in paper III.) 

The outline of this paper is as follows: In section y 
we introduce the NTKM, and in section III outline our 



culated in section ^ Our results for T{v) are compared 
to experiment and the NCA results of Hettler, Kroha and 
Hershfield in s ection [V]|, and our conclusions summarized 
in section VII. 



strategy for solving it by a combination of CFT methods 
with Hershfield's F-operator approach. This strategy is 
implemented in section IV, where the scattering states 



are calculated. The current and scaling function are cal- 



More than half of the paper is taken up by appen- 
dices. The lengthier ones (^qJqQ^,^) summarize, for 
the sake of convenience, background material that is as- 
sumed known in the main text; the others (|^,|G|J|) contain 
original work related to the main text. In Appendix H, 
we recall some standard results from the semi-classical 
theory of non-equilibrium transport through a ballistic 
nanoconstriction. Appendices H and provide a brief 
review of the recent series of papers by Zarand and Za- 
wadowski on the (bulk) non-magnetic Kondo model and 
its renormalization toward the 2CK model at low tem- 
peratures. Recent criticism of their conclusions are dis- 
cussed in Appendix y. In Appendix ^ we compare our 
CFT results with those from the poor man's scaling ap- 
proach in the limit of large channel number (fc — *■ cxd), 
in which the latter approach becomes exact. Hershfield's 
y-operator formalism is briefly reviewed in Appendix ^. 
Appendix O illustrates the general scattering states for- 
malism developed in sections IV and M with an simple 
example. In Appendix H we give some background on 
the NCA calculations of Hettler, Kroha and Hershfield. 
Finally, in Appendix | we discuss F/T^-correction to our 
results. 



II. THE NANOCONSTRICTION 
NON-MAGNETIC KONDO MODEL 

In this section we introduce a new model, to be 
called the nanoconstriction two- channel Kondo model 
(NTKM), to describe the interaction of conduction elec- 
trons with a TLS in the nanoconstriction. We shall take 
as guideline the results of Zawadowski and coworkers, 
who introduced the non-magnetic Kondo Hamiltonian 
to describe the TLS-electron interaction (summarized in 
Appendix H) and showed that under renormalization it 
flows towards the non- Fermi- liquid fixed point of the 2CK 
model (in a way summarized in Appendix 0).^ However, 
we shall not be interested in the details of the renormal- 
ization process from some bare to some effective model. 
Instead our attitude here is that of phenomenologists: 
since the detailed microscopic nature of the presumed 



^It should be pointed out that the question as to 
whether a realistic TLS-electron system will reach the 2CK 
non- Fermi- li tjuid .jesipie under renormalization is currently 
controversialoOEj'Ed (see Appendix pi). In the present pa- 
per, though, we do not attempt to clarify any of the con- 
troversial issues. We simply take the view that it would be 
useful to know what the scaling curve would look like if the 
system indeed does reach the 2CK non-Fermi-liquid regime, 
and hence do the calculation, assuming it does. 



TLSs is unknown, so too is the "correct" microscopic, 
bare Haniiltonian. The best one can hope for is to find 
a phenomenological Hamiltonian that satisfactorily ac- 
counts for the observed phenomena. As argued at length 
in paper I, the 2CK model with energy splitting A ~ 
passes this test on a qualitative level. We regard this 
as sufhcient justification to use 2CK ideas as a basis for 
quantitative calculations, in order to test whether quan- 
titiative agreement with experiment can be achieved. 

The NTKM that we shall write down is the simplest 
model we can think of that contains the non- Fermi-liquid 
physics of the 2CK model, but also accounts for the com- 
plications brought about by a nanoconstriction geometry 
relative to the bulk situation. We introduce it as a phe- 
nomenological Ansatz, without attempting to provide a 
detailed microscopic derivation. Since our aim is to cal- 
culate a universal curve, characteristic of the 2CK model 
but experimentally found to be sample-independent, we 
believe that such lack of attention to microscopic details 
has experimental justification. 

The main complications arising in a nanoconstriction 
geometry relative to the bulk case are, firstly, that one 
has to distinguish between electrons leaving and entering 
the L and R leads, and secondly, that the application of 
a voltage induces a non-equilibrium electron distribution 
in the nanoconstriction. 

We thus have to deal with a non-equilibrium problem 
with non-trivial interactions,-. The standard procedure 
(due to Kadanoff and BaymE3) for defining such a prob- 
lem requires conceptual care and may for clarity be or- 
ganized into six steps: 

First the problem is defined in the absence of inter- 
actions, by defining 

(51) a free Hamiltonian Ho with eigenstates {|£?7)o}, 

(52) a free density matrix po governing their non- 
equilibrium occupation, 

(53) and the physical quantities of interest, in our case 
the current / (with expectation value 

(/) = Tipol/Trpo in the absence of interactions). 
Then the interactions are switched on, by defining 

(54) an interaction Hamiltonian Hint, 

(55) and the full density matrix p, which governs the 
non-equilibrium occupation of states 



for the fully interacting system. (Typically, this is done 

by adiabatically switching on 

Hint, and keeping track of how the initial po develops 

into a final p.) 

(S6) Expectation values are calculated according to (/) = 

Trp//Trp . 

In this section, we address steps (SI) to (S4). [(S5) 
and (S6) are discussed in sections IV ^, respectively]. 
We also explain, within the poor man's scaling approach, 
why the flow towards the non-Fermi-liquid regime is not 
disrupted by y 7^ as long as F <C T^- 



A. Free Hamiltonian Ho 

We consider a single TLS at the center of the nanocon- 
striction (see Fig. 1 of paper I for a scetch of the nanocon- 
strictions used in the RB experiment). We consider only 
those modes of electrons that contribute to the ZBA, i.e. 
that interact with this TLS when passing through the 
nanoconstricion. 

To describe these electrons, we imagine that the "free 
nanoconstriction Schrodinger equation" for free electrons 
and some random static impurities but no TLS-electron 
interaction, with boundary conditions that all electron 
wave- functions vanish on the metal-insulator boundary, 
has already been solved (impossible in practice, but not 
in principle). This provides us [step (SI)] with a com- 
plete set of single-particle eigenstates {\e,ri)o = cj,g^|0)} 
(where |0) = vacuum), in terms of which Ho is diagonal: 



Ho 



V 



D 



de e cL 



(2) 



Here the continuous energy label e is taken to lie in 
a band of width 2D, symmetric about the equilibrium 
Fermi energy (at e = 0) , with constant^ density of states'^ 
No- The latter has been absorbed into the normalization 



of the cjjg^'s, which we take as 



{c, 



06171 (^oe-qi 



5r,r,'S{e - e') 



(3) 



^Very recent work by Zarand and UdvarditJ has shown that 
using a constant density of states is probably less realistic in 
a nanoconstriction than in the bulk (where it is standard), 
because the local density of states fluctuates strongly as a 
function of r and e. This is the kind of complication that our 
phenomenological approach has to ignore. 

^Since the density of states diverges for inflnite systems, 
the expectation values of some operators, e.g. the current 
[e.g. see footnote 5 and Eq. (B9|)], have to be evaluated in 
a finite system with a discrete energy spectrum. In such 
cases, we use the replacement rules: Jde — > N~^ ^ , and 
5{e - e') -^ NoS,,, . 



The label rj collectively denotes a set of discrete 
quantum numbers, rj = {a, a, i) = (species, pseudo- 
spin, channel)-index, which have the following meaning: 
i ^tii is the electron's Pauli spin, which will be seen 
below to play the role of channel-index in the NTKM. 
a = 1,2 is a discrete pseudo-spin index, the nanocon- 
striction analogue of Vladar and Zawadowksi'siiangular" 
index a [see e.g. Eq. (2.36) of the first paper ofo; irill3, a 
was called a "parity" index]. It labels those two sets of 
free states {|e, cr, l,i)} and {|e,cr, 2, i)} that in the non- 
Fermi-liquid regime will couple most strongly to the TLS. 
For example, if the free wave- functions were expanded in 
terms of angular harmonics, a = 1,2 would label two 
complicated linear combinations of y;_m(0, (/>) functions. 
Strictly speaking a can take on a large number of discrete 
values, but we ignnte all but two, in the spirit of Zawad- 
owski's bulk resultllj that the others decouple from the 
impurity when the temperature is lowered and the sys- 
tem flows toward a non- Fermi-liquid fixed point with an 
effective electron pseudo-spin of i. (The modes we ig- 
nore contribute to the background conductance, but not 
to the ZBA.) 

Finally, a = (+, — ) = {L,R), the species index, de- 
notes the direction of propagation of the incident elec- 
tron: cr = L = + for left-moving electrons, incident to- 
ward the left from z = -|-oo in the right lead; a = R ^ — 
for right-moving electrons, incident toward the right from 
z = — oo in the left lead. (For example, in spherical 
coordinates the asymptotic behavior of the incident (or 
transmitted) parts of the wave-function of both L- and 
_R-movers will be proportional to e~''^''/r (or e^^"^ /r) as 
r — > cx).) The nanoconstriction geometry necessitates 
this distinction between L- and i?-movers (not needed in 
the bulk case) , firstly because L- and _R- movers originate 
from different leads, which are at different chemical po- 
tentials if V^ 7^ 0, and secondly because they contribute 
with different sign to the current. 



^) of -\-eV/2 and —eV/2, respectively.^ As input, we 
use a standard result from the semi-classical theory of 
non-equilibrium transport of electrons through a ballis- 
tic nanoconstrictioiiE3 (summarized in Appendix |A|) : At 
the center of the constriction, the distribution of occupied 
electron states in momentum space is highly anisotropic 
(see Fig. in Appendix H). It consists of two sectors, 
to be denoted by L or R, that contain the momenta of 
all electrons that are incident as L or R-niovers, i.e. are 
injected from the R oi L leads. Consequently, the Fermi 
energies of the L/R sectors are equal to those of the R/L 
leads, namely /i± = zL^eV. 

We formalize these standard results by associating the 
L/R sectors with the species quantum number a = 
L/R = ± introduced above (correspondingly /i,j will 
stand for fJ-±), and adopting the following form for the 
free density matrix po'- 



p^^,~m.-y.\ 



(O)o ^ ^ , (4) 



Trpo 



where the Yo-operator is defined by 



Yo = \eV (Nl -Nr) ^Y. i"^ J'^^ 4„Co 



(5) 



Here Nl and Nr denote the total number of L- and R- 
moving electrons.^ It follows that (in the absence of in- 
teractions) 

(4„(r)c„,,,,(r')) = e'^^-^'^f{e,r,)S^^,S{e~e') , (6) 
where /(e,^) ^ -^^^-1^^ . 



C. The free current through the nanoconstriction 



B. The free density matrix po 

We now turn to step (S2), the definition of po, the 
free density matrix for Hint = but arbitrary volt- 
age. The right and left leads have chemical potentials 
(measured relative to the equilibrium chemical potential 



The ZBA arises from backscattering by the TLS of 
electrons that would otherwise have passed through the 
constriction. Thus, we assume that they would con- 
tribute one unit e^/ft. of conductance if the interaction 
were turned off. (More generally, one could use 7^e^//i, 
where 7^ is a transmission coefficient, but this only af- 
fects the (non-universal) amplitude of the ZBA.) Thus, 



''Our figures and arguments are given for the case eV > 0. 
We take e = -|e| and hence V = -\V\. With ii± = ±eV/2 
for R/L leads, there then is a net flow of electrons from right 
to left, and the current to the right is positive. 

^To evaluate (cJejjCo^/^) for e = e' , we have to give mean- 
ing to S{e — s) of Eq. (h), which seems to diverge because 
we took the thermodynamic limit of an infinitely large sys- 
tem. We do this by replacing S{e — e) by the corresponding 
finite-system expression of NoS^gri [see footnote 3], i.e. we use 



\<--0£ 



,) = f{e,v)No 



we may define [step (S3)] our current operator simply as 
the difference between tire number of electrons transmit- 
ted as L- or _R-movers: 



/ = 






deer cj,.„c 



oeri'^oeri 



(7) 



Our signs are chosen such that (/)o > if the net flow of 
electrons is from right to left, while the prefactor \e\/hNo 
is needed, because of our choice of normalization, to 
obtain^ a conductance of e^/h per channel. 



studied by Hettler et. al&B&B using numerical NCA tech- 
niques, with whose results we shall compare our own (see 
section 



VIC) 



The Hamiltonian introduced above is strictly speak- 
ing not a 2CK Hamiltonian, since a = zt and i =1, | give 
four different combinations of indices that do not Kondo- 
couple to the impurity. However, it can be mapped onto 
a 2-channel model by making a unitary transformation. 



^oer] — ^^ rjrj^oerj 7 



Nf^rj 



-^^ (7(7^aa.^ii 



(10) 



D. The nanoconstriction 2CK interaction 



chosen such that it diagonalizes v^ 
choice (||) for Vaa' , N^a- is given by 



For our present 



We now come to step (S4), the specification of the 
electron-TLS interaction, for which we make the follow- 
ing phenomenological Ansatz: 






Vr 



rjrj' — '^K'^crcr'^ii 



' Oii' \ 2^aQ. 



rjrj' Coe'7]' ; 



(8) 



A^^ 



^(i-i)„^ {»"»-%. = (; s)„, 



V2 



(11) 



We shall refer to the operators Coerj as L/R operators and 
the Coef] as even/odd operators, and always put a bar over 
all indices and matrices refering to the even/odd basis. 
In the even-odd basis, the interaction becomes 



Here S is the TLS pseudo-spin operator acting in the 
two-dimensional Hilbert space of the TLS. Following the 
assumption (A2) of section V.D of I, we henceforth as- 
sume that A, the TLS excitation energy, is the smallest 
energy scale in the problem, and set A = 0. 

As far as the pseudospin and channel indices a and i 
are concerned. Hint is simply the isotropic 2CK Hamilto- 
nian to which, according to Zawadowski's analysis for a 
bulk system, a realistic TLS coupled to electrons will flow 
at sufficiently low temperatures. However, we introduced 
an extra Hermitian 2x2 matrix Vaa' , which enables an in- 
cident electron, say a L-mover, to be scattered into either 
a L- or a i?-mover, independent of whether its pseudo- 
spin index a and that of the TLS do or do not flip.^ In 
general, v^^' can be any Hermitian matrix, but, for rea- 
sons given below, it is actually sufficient to consider only 
the very simple case 



*(1 



(9) 



Note that with this choice, our model is equivalent (af- 
ter a Schrieffer-Wolf transformation) to a model recently 



Hint = de /rfe'^ cl^ijCoeTf V 



(12) 



V nri' — VkOhi 







Thus, in the even/odd basis, one set of channels, the 
odd channels {a = o), completely decouples from the 
impurity. The other set of channels, the even channels 
(ct = e), constitute a true 2CK problem, which will even- 
tually be responsible for the non-Fermi liquid behavior 
of the NTKM. 

If one chooses a more general v^a' than Eq. (S), the odd 
channel will not completely decouple, but (barring some 
accidental degeneracies) the even and odd channels will 
always couple to the TLS with different strenghts. At low 
enough temperatures, the one coupled more weakly cm 
be assumed to d ecouple completely (a la Zawadowskill3, 
see section |C3| of Appendix g), leaving again a 2CK 
problem for the even channel. This is the reason why it 
is sufficient to take Vca' as in (||). 



® Note that the interaction of Eqs. (H) and ffih is reminiscent 
of the tunneling Hamiltonian Htun in the standard problem 
of electrons tunneling through an insulating barrier that sep- 
arates two electronic baths: the oflt-diagonal components of 
Vera' transfers an electron from one bath to the other, with 
the implicit assumption that this does not disturb the thermal 
distribution of electrons in the baths significantly. 



E. Poor Man's Scaling Equations unaffected by V 

The model we wrote down assumes that the NFL 
regime of the TLS-elcctron system has aheady been 
reached. However, one may wonder whether having 
y ^ would not prevent the TLS-electron system from 
reaching the non-Fermi-liquid regime at all. That this is 
not the case for V sufficiently small (<C Tj^) can be seen 
by the following poor man's scaling argument: Since the 
poor man's scaling equations are derived by adjusting the 
cut-off from D to D' , which are both ^ V,T, they are 
independent of V for the same reason as that they are 
independent of T (namely the change in coupling con- 
stants needed to compensate D ^ D' does not depend 
on energies V and T that are much smaller than D). 
In other words, the scaling equations for V ^ are the 
same as those for V = 0, meaning that the initial RG 
flow is unaffected by F 7^ 0. Eventually, the RG flow 
is cut off by either V or T, whichever is larger, result- 
ing in an effective Hamiltonian that depends on V or T. 
However, if both are <C T^, the RG flow will terminate 
in the close vicinity of the non- Fermi- liquid fixed point, 
even li V 7^ 0, and the V or T-dependence of the effec- 
tive Hamiltonian will be of order V/T^ or T/T^, both 
<^ 1. This is the basis of our key assumption, stated 
in the introduction and implicit in the Ansatz (g), that 
for V" ^ but V/Tk <C 1, the non- Fermi- liquid regime 
is governed by essentially the same effective Hamiltonian 
as for V = i). 



III. OUTLINE OF GENERAL STRATEGY 

We now have to address step (S5) of the process of 
defining a fully interacting, non-equilibrium problem, 
namely the definition of the full density matrix p for 
V ^Q and Hint 7^ 0. In this section, the heart of this pa- 
per, we propose a strategy for doing this which combines 
ideas from CFT with Hershfield's F-operator formulation 
of non-equilibrium problems. The section is conceptual 
in nature; technical details follow in sections IV and M, 
and in paper HI. 



A. Hershfield's Y-operator approach to 
Non-Equilibrium Problems 

Typically, the full p is defined by adiabatically turning 
on Hint and following the evolution of the initial den- 
sity matrix po to a final p (see Appendix ^) . Expanding 
the time-evolution operator in powers of Hint , one then 
generates a perturbation expansion that can be handled 
using the Keldysh technique. 

However, for the Kondo problem, perturbation the- 
ory breaks down for T < T^, where many-body ef- 
fects become important. Therefore we shall adopt 
Hershfield's so-called [^-operator formulation of non- 
equilibrium problemsE-3, which is in principle non- 
perturbative. 

The main idea of Hershfield's approach (briefly sum- 
marized in Appendix ^, is as follows. As the interaction 
Hint is adiabatically turned on, the density operator adi- 

-/3(H„-n) 

~0(H-Y) 



abatically evolves from its initial form po — 
into a final form that Hershfield writes as p = e 



This defines the operator Y , which is the adiabatically 
evolved version of Yq and is conserved ([F, H] =0). The 
formal similarity between p and po implies that when ex- 
pressed in terms of H and Y , the non-equilibrium prob- 
lem has been cast in a form that is formally equivalent 
to an equilibrium problem. 

This becomes particularly evident if one considers the 
set of simultaneous eigenstates of H and Y , which we 
shall call the scattering states and denote by {\eri) = 
c\^\Q)}. Loosely speaking, they can be viewed as the 
states into which the free basis states {|e?7)o} develop as 
Hint is turned on (in the sense that c\n is some function 

of the {Coe/„/}, which reduces to cj^^^ for Hint = 0). For 
scattering problems like the NTKM, in which a free elec- 
tron is incident upon a scatterer and scatters into some- 
thing complicated, there evidently must be a one-to-one 
correspondence between the states \eTf)o and \eri): the 
incident parts of their wave- functions {x\eri)o and {x\erf) 
must be identical. (The outgoing parts, which contain 
scattering information, will of course be different - this 
will be made explicit in Eq. ( p4| ) below.) This is why the 
free and scattering states can be labelled by the same 
indices, and also have the same density of states.^ 

Furthermore, for such scattering problems, H and Y 



'^One might ask whether the very notion of scattering states 
make sense for a dynamical impurity problem, since the scat- 
terer is constantly flipping its pseudo-spin. However, in the 
CFT solution of Kondo problems, the impurity completely 
disappears from the scene (being absorbed in the definition 
of a new spin current, see paper III) . Thus the theory contains 
only electron degrees of freedom, for which one can meaning- 
fully introduce scattering states. 

*For problems other than scattering problems, Eq. (rMj) does 
not necessarily hold. 



will have the following form:' 



H = Y^ descl^Cer, , (13) 

V •' 

Y = Y, [det^^cl^c,^ i^Yo). (14) 

The form used here for Y follows because Y evolves 
from Yo as Hint is turned on, implying that Y can be 



obtained from Yo by replacing the Coer, in Eq. 



by 



the scattering-state operators Ce,, into which the latter 
evolvetZl. Eq. (03) and mm imply that non-equilibrium 
thermal expectation values of the c^^'s have the standard 
form: 

(4,(r)ce,,,(T')) - ef^^-^')f{e,^)5,„.5{e - e') (15) 

where fis,"!]) 



p/3(£-A'.,) + 1 

This is precisely the same form as that satisfied by the 
non-interacting Coer/'s in the absence of interactions [see 
Eq. (|6|)]. The intuitive reason for this remarkably simple 
result is clear: the Boltzman weight of a scattering state 
must be the same as that of the corresponding free state, 
since the thermal equilibration that leads to the Boltz- 
mann factors happens deep inside the leads, before the 
electrons are injected and scattered by Hint (this of course 
remains true when L- and R leads have different chemical 
potentials - all that happens for y 7^ is that the oc- 
cupation probabilities pick up a ^-dependence reflecting 
from which lead the electron was injected). 

This result provides us with a very clear picture of how 
the current through a nanoconstriction should be calcu- 
lated: when injecting electrons from the leads into the 



constriction, the thermal weighting is done precisely as 
for free particles, i.e. an electron incident in the state 
\e'rj')o is injected with weight f{e',rj'). For each such 
electron, one has to determine the scattering amplitude 
Ur/r/'is'), i-e. the amplitude with which it emerges from 
the scattering process in the state \s',ri)o (where we as- 
sumed elastic scattering). These ampli tudes (defined 
more explicitly below, see section [V C ) are the non- 
trivial ingredients of the scattering states, which contain 
all relevant information about the scattering process.^ 
Once they are known, it is straightforward to calculate 
the current as a thermally weighted sum over transmis- 
sion probabilities. 

Since expectation values expressed in terms of scatter- 
ing states are so simple, it is useful to reexpress all phys- 
ical operators in terms of them. To this end, we define 
Uj-i'nie' , e) = ois'v'l^v) to be the unitary transformation^ 
that relates the scattering states to the free basis states: 



(16) 



k'?)^E/^^'l-^V)oC/,',(e',£) 

(5„„.,5(£ - e') = E \^^ f^^y^' ^~)^'>'' (^^' ^') ■ (18) 

For example, the current of Eq. (0) takes the form: 



No 



xRe 



[aUl,^{e',e)U,^.{e,e"){cl,^,Ce.^. 



For a many-body problem such as the Kondo problem, 
complicated combinations of particle-hole excitations are cre- 
ated upon scattering, which can not simply be written as a 
linear combination ^ c^_, ?7^^/(e') of single-particle excita;- 
tions. Hewever, it was shown by Maldacena and Ludwi^3 
(see alsocfl) that the scattering matrix for free electrons in- 
cident on a Kondo impurity is unitary if the single-particle 
Hilbert space of free-electron states {|£?7)o} is appropriately 
enlarged to include "Kondo excitations" (called "spinors" by 
them, see paper III). This means that the outgoing states 
can be written as linear combinations of free-electron states 
{k'?)o} and a new set of Kondo excitation states {|e??)s}. 
The corresponding set of creation operators {co^^} are com- 
plicated, highly non-linear functions (not mere linear com- 
binations) of the {c^ , } and will be discussed in paper III. 
Thus, in the formalism developed below, the unitary trans- 
formation in Eq. (Ilq) is implicitly understood to act in the 
enlarged Hilbert space of {\eTf)o, |£'7)o} states, and the collec- 
tive index r] implicitly includes another index o = (/, k) to 
distinguish free from Kondo states. However, this will only 
be made explicit in paper III. 



The reality of / is of course automatically ensured by the 
hermiticity of the current operator, and the reminder Re[ 
] has been inserted merely for future convenience. 

We shall show below that the [/,,,,'(£,£'), and hence 
also the current, are completely determined by the 
Ujjri'is')- Unfortunately, Hershfield's formalism gives no 
recipe for finding these explicitly for a given problem. 
Thus, the crucial question now becomes: how does one 
calculate the scattering amplitudes? 



B. Equating CFT- and scattering-state Green's 
Functions 



In general, finding the scattering amplitudes is just as 
difhcult as solving the problem by other (e.g.. Keldysh) 
methods. (For example, Keldysh methods are used to 
find the scattering states ofja-closely related Kondo prob- 
lem in Appendix C of Ref.li3.) However, for V = the 
even sector of the NTKM is equivalent ta.thp 2CK model, 
which AL solved exactly using CFTBElkrEl. (This 
equivalence is shown explicitly below, when we rewrite 
the model in field theoretical language, see Eq. ( p7|) and 
( PSJ ) below.) Therefore, we propose that the scattering 
states of the NTKM can be extracted from AL's results. 
We now explain how this can be done. 

One of AL's central results is an explicit and exact 
expression for the equilibrium Green's G^^' = ^{'i'TivL') 
[defined explicitly in Eq. (p9|)], which gives the ampli- 
tude that an incident ?7'-electron will emerge from the 
scattering process as outgoing //-electron. Evidently, it 
must contain information about the scattering ampli- 
tudes. Indeed, we shall show that when the same equilib- 
rium Green's function is calculated explicitly using the 
scattering state formalism, it is completely determined 
by C/^^'(e'). Therefore, by equating the scattering- states 
form for Gr^rf' to the corresponding CFT result, C/^^/ (e') 
can be extracted from the latter. 

Of course, this procedure only yields the V = Q value 
of Urirt' , whereas to calculate the nonequilibrium current, 
we actually need its F 9^ values too. Moreover, it is 
clear that in general C/^^' must depend on V, since if V 
is sufficiently large, it is known to non-trivially affect the 
many-body physics of the Kondo problem. For example, 
for V j^ 0, the difference in Fermi energies of the L- and 
R leads-causes the Kondo peak in the density of states to 
splitESEj into two separate peaks (at CTiergies /i ± ^eV, 
see Fig. H of Appendix H, taken fronJlj). Moreover, the 
effective Hamiltonian in poor-man's scaling approaches 
depends on V if it is th e lar gest low-energy cut-off in 
the problem (see section HE), and if V is too large, it 
will cut off the renormalization group flow towards the 
non-Fermi-liquid fixed point before the non-Fcrmi-liquid 
regime is reached. 

However, such ^-induced effects should be negligible 
for sufficiently small V. For example, when V -^Tk, the 
splitting of the Kondo peak by eV is negligible compared 



to its width, which is ex T^. Said in poor-man's scal- 
ing language, if (T <) 1/ < T^^ , then V ^0 cuts off the 
renormalization group flow at a point sufficiently close to 
the non-Fermi-liquid fixed point that the physics should 
still governed by the latter. Hence, wc propose that in the 
non- Fermi- liquid regime oi V <ti T^, the T^-dependence 
of the scattering amplitudes C/,,^' is negligible, and hence 
shall always use their V = values below. (In a sense, the 
condition that this procedure be valid can be regarded as 
our definition of the "non- Fermi-liquid regime".) More 
formally, we assume that f/^^' can be expanded in powers 
of V/Tfc , and use only the zeroth term. (In Appendix |, 
we show that the leading V/Tk correction only produces 
a subleading correction to the desired scaling function.) 

The intuitive motivation for neglecting the In- 
dependence of the scattering amplitudes is based on the 
assumption that the effect oiV^O can be characterized 
as follows ii V -^ T^: although the leads inject electrons 
into the non- Fermi-liquid state that, since T^ ^ 0, are 
able to probe its nature at energies different from Ef, 
they only probe gently, i.e. they inject sufficiently few 
that the non-Fermi-liquid state itself is not disrupted. 
Since the "output" of this probing, namely the scatter- 
ing amplitudes, depend non-linearly on e, the current will 
depend non-linearly on V, too, even if t/^^/ (e') itself is 
F-independent. 

Another underlying assumption of our proposed strat- 
egy is that the strong-coupling or fixed-point fields 
ipri (t, ix) occuring in the CFT treatment can be expanded 
in terms of a set of fermionic excitations (though these 
are very complicated non-linear combinations of the free 
ones, cf. footnote 9), else it would not make sense to 
equate a CFT Green's function to one constructed from 
scattering states. That this is indeed the case will be 
shown in paper III. 



IV. EXTRACTING SCATTERING STATES FROM 
CFT RESULTS 



To implement our strategy for finding Ur/n', the first 
step is to rewrite the NTKM of section O in field theory 
language by introducing a set of fields ipri (ix) . Then we 



define the Green's function G„ 



<T^P,^Pl, 



and show 



that it is completely determined by Ur/n' (which turns out 
to be its spectral function). Finally, we equate this Gr/n' 
to the corresponding exact CFT result of AL, which al- 
lows us to obtain the corresponding exact expression for 
Urjri' explicitly. 



A. Transcription to Field Theory 

To rewrite the "bare" NTKM introduced in section || 
in field theory language, we introduce for each channel 



r] a 1-dimensional, second-quantized field ijJniTjix) (with 
X e [— /, I], I ^ oo) as a Fourier-integral over all e:^° 



■ipniix) 



d£ e-'^^/'"'^ c. 



osrj ; 



-oer] 



1/2 



1 K^ / ||e'-/'^-^,M, 



;/2 



{'(/'j,(«a;),'0l'(ia::)} ^2'k6,^,^'5{x ~ x) 



(20) 

(21) 
(22) 



The factors of 1% and f p, inserted for dimensional reasons, 
are henceforth set = 1. 

Note that "ipriiix) is not the usual electron field '^{x), 
which is constructed from the actual (unknown) wave- 
functions {x\er])o through ^(x) = '^ jde{x\eri)oCoeri- In- 
stead, ipriiix) is best thought of simply as the Fourier 
transform of Coeri , this being a convenient way of rewrit- 
ing the problem in field-theoretical language. Neverthe- 
less, the role of x is strongly analogous to that of the 
"radial" coordinate of the actual wave- function \I>e^(af), 
and ipli (x) can be interpreted as the operator that creates 
an electron with quantum numbers 77 at ''position" x. 

Using Eq. (|2l|), Ho and Hint of Eqs. (||) and (||) can 
be written as 



/oo 

77 jj' 



(23) 
(24) 



By simple Fourier transformation, we have hence ar- 
rived at a 1+1-dimensional field theory, defined by 
Eqs. ( p3| ) and (p5). The reason why this (and not a 3-1-1 
dimensional theory) resulted, is essentially that there is 
only one continuous quantum number, namely e, in the 
problem, with respect to which we can Fourier trans- 
form. This in turn is a result of the constriction ge- 
ometry, which defines a definite and unique origin, and 
consequently a notion of a single "radial" coordinate (in 
spherical coordinates it is the radius r), to which our 
X roughly corresponds. Moreover, the fact that we as- 
sumed a constant density of states and hence a linear 
dispersion implies that the free fields are conformally in- 
variant, which is the key property required for the sub- 
sequent application of AL's CFT methods. 

The Heisenberg equation of motion. 



dril^nir, ix) = [i/'.,(t, ix),Ho + Hint] 

= {S^pyidx + 2TT5[x)V^^>)ipr)'{ix) 



shows that for all a; ^ 0, the fields depend only on r + ix. 
[This is the reason for writing the argument of V'r; as (ix) 
in Eq. (|2^) , since the r dependence of tpn can then simply 
be obtained by analytic continuation (ix —^ r+ix).] Con- 
sequently, by construction, all fields are "mathematical 
left-movers" , incident from x = cxd and traveling toward 
X — — cxD. The effect of the scattering term Hint is to mix 
the different incident channels with each other at a; = 0, 
so that tpri{T, ix) will differ from a free field only for a; < 0. 
Thus, we have turned our problem into a one- dimensional 
scattering problem, with all free fields incident from the 
right, and all scattered ones outgoing to the left. This is in 
exact analogy to AL's treatment of the Kondo problem, 
which in fact was the motivation for introducing both 
physical L- and i?-movers as "mathematical left-movers" 
in Eq. (20). Of course, the distinction between physical 
L- and i?- movers is carried by the index a — L, R, and 
L-R backscattering is described by the a ^ a' terms in 

V rin' . 



B. Transformation to even-odd basis 



As mentioned in section II D, the relation between the 



NTKM and the standard 2CK model is best understood 
in the even-odd basis (denoted by bars) of operators 
Coefj = NfjrjCoerj [scc Eq. (|l^]. Therefore, we define even- 



odd fields 



i^fiiix) = Nf^r)i^n(ix) , 



(26) 



(25) 



normalized as in Eq. (p2|). In terms of these. Ho and Hint 
of Eqs. (H) and (||) are: 

/°° _+ _ 

!^i?^{zx)id^i,^{zx) , (27) 

,1 -°° 

H„u = Y. ^^ai(0) (y^5,,\cj^^, ■ S) ^^a^(0) . (28) 



The odd channel {a ~ o) decouples from Hint- In the 
even channel {a — e), Ho + Hint is precisely the "bare" 
Hamiltonian of the ecuiilibrium 2CK model solved ex- 
actly by AL [see e.g£3, Eq. (2.17)]. Therefore, the 
even channels will display non-Fermi-liquid behavior for 
T,V -^T^. 



^"strictly speaking, the Jde integrals have to be cut oflt, 

nD 

j_ de, at a bandwidth D satisfying T,V <^ D. However, we 
take D —* 00 (since the errors thus introduced are of order 
T/D, V/D <C 1 and hence negligible even for finite D). This 
allows us to invert relations such as (GCl) straightforwardly. 
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C. Definition of scattering amplitude [7^^/(e') 

Having rewritten the model in field theory language, 
we can define the equilibrium Green's function (in the 
original L-R basis) that is to be the link to AL's CFT 
results: ^^ 

G^^,iT,-ir;T',ir') = -(rV^(r, -zr)Vj,(T', ir')) , (29) 

with r,r' > 0. Since its arguments correspond to taking 
X = — r < and x' = r' > 0, it gives the amplitude that 
an incident ?7'-electron will emerge from the scattering 
process as outgoing yy-electron. 

In order to calculate Gr^r]' in terms of scattering states, 
we rewrite the fields ■0?) (t, ia;) in terms of the c^^'s. In- 
serting the inverse of Eq. (n7h into Eq. (20) and defining 



<pe'v'{ix,ri) = de 



se-'"U^,yie,e') 



we find 



1p^{T,ix) = > lde'(J3e'n'{ix,ri)Ce'r,'(j) 



(30) 



(31) 



Now, because scattering takes place only at a; = 0, for 
X > (i.e. before the scatterer is encountered) the wave- 
function 0£'^'(ia;, 7y) must correspond to the free wave- 
function e"'"^ ^ of the state \£'rj')o- Thus, we make the 
following Ansatz:^'^ 



(l3^>ri'{ix,ri) 



^~ie X 



U^^,ie')e{-x)+5^^,e{x) . (34) 



This relation defines the matrix C/^^/(e'), which clearly 
can be interpreted as a scattering amplitude, since it 
specifies the amplitude for an electron incident with 
quantum numbers {e'rj') to emerge with quantum num- 
bers {e'r]). 

The relation between the scattering amplitude Ur,r]' (&') 
and the matrix {/^^/ (e, s') can be found by inserting 
Eq. (H) into the inverse of Eq. (|^): 



t/,,. (£,£')= /|fe--(/.e','(za;,r;) 



1 

27ri 



e — e' — it 



^nv' 



e- e' 



(35) 
(36) 



which implies that <j)^',^'{ix,ri) = {ipri{T,ix)cl, ,{t)). 
Since by its definition (pG) ■0jj(r, ix) has the interpre- 
tation of creating an electron with quantum numbers i] 
at x, this shows that ^^'^/(ia;, 77) may be thought of as 
the "wave- function" for the scattering states je'ry'):^^ it 
gives the amplitude for an electron in state \e'r]') to be 
found at x with quantum number 77. The orthonormality 
and completeness of these wave-functions is guaranteed 
by the unitarity (|l|) of f/^,,'(£, e'): 



H / If 't^lr,{ix,fi)<l)e'y'{ix,f}) ^5^^'5{e~e' 



y^ Ide 4>lfi{ix,ri)(j)if,{ix' ,r]') = 27r (5,,,,/(5(a; - x') 



(32) 
(33) 



(e > is infinitessimally small). This shows that 
Uri,f{s,e') is completely known once [/,,,,'(£') is known. 
The unitarity condition Eq. dlj) on C/^^/ (e, e') then im- 
mediately implies unitarity for U^py (e') (the Jde integral 
can trivially be done by contour methods): 






(37) 



The unitarity of U^fi{e^ could of course also have been 
anticipated from Eq. (p4|): it ensures that scattering con- 
serves probability, i.e. that ^ |(/<j'^'(ia:, 77)^ is the same 
for a; > and x < 0. 

The current can be rewritten as follows by inserting 
Eq. (H) into Eq. (|l|): 



rjrj' T]" 



^ y Ide' I d£"Re 




e' - e" - 2ie 



0^'^0'qri" 



J 



e' - £'• 



2ie 



(c\ , 



(38) 



In paper III, this Green's function is denoted by 
G^'^/{z*;z'), following the notation used AL. 

^^This interpretation of (j)^in'{ix,ri) as a wave-function is 
meant as a mnemonic and should not be taken literally; as 
mentioned in section II A, the actual physical wave-functions 
are intractably complicated. 

^^In writing Eq. (34), we have assumed elastic scattering 
{sin — Eout)- For a 2CK model, this holds only if the impu- 
rity energy splitting A = 0, as assumed in this paper, so that 
electrons cannot exchange energy with the impurity. 
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Ul4e')U^^4e") + ^,',V' 5{^ - n (cl'^'Ce",") 



^ /deVi \ul,^{e')U^^'{e') + 5., J f{e',v') 



~ h 

n'v 

To obtain Eq. (|3|), the Jde integral in Eq. (|l|) 
was done using contour methods. Eq. ( p^ ) fol- 
lows since the diagonal nature of {cli,,„,Coe"ri") en- 
sures that UL {e')Urj-q"{e") is real, so that we may use 

Re [(27ri)(£' - e" T 2ie)]^^ = ±iJ(e' - e"). Finally, to 
obtain Eq. (EG), we used footnote 5. The problem of 



(39) 
(40) 



finding the scattering amplitude C/^^/ (e) . 



D. Extracting [7,,^/ from the Green's Function G,, 



calculating the current has thus been reduced to that of Eq. (g9|) can be reduced to the form 

G^^>{t, -ir; t', ir') = - Y. j de de' U^f,{e)6f,,^, (c,^(r)4,-,(T')) e-'^-'^-''"'"' 



Using Eqs. (|l|), (||) and (|l|) (with ^l' ^ 0), G^^. of 

(41) 



de U^yy (e) 



_, — e(r— -ir— r' — 'ir') 

e-P^ + 1 



Its Matsubara-transform is readily found to be 



G,,,f{iujn;r,r') = de 



_t^^^.(e)e-('-+'-') 



lOJn 



(42) 



(43) 



This is the central result of this section: G^pf is completely determined by U^,j' (e) , which is proportional to the spectral 
function of G,,,/, as is evident from the form of Eq. (H3). Conversely, if Griri' is known, C/^^'(e) can be extracted from 
it using 



Urin'ie) 



_ _j_g-ie(r+r') 

27T 



[G^^,{e + iO+-r,r') - G,„,{e - iO+;r,r')] 



(44) 



Thus, by equating Gr/n' to the corresponding exact CFT in paper III. For present purposes, it suffices to consider 
result, if,„'{e) can be extracted from the latter using CFT as a "black box" that, starting from Eqs. (|^ and 
Eq. (0;. In the next section section, we shall cite the (H), allows one to calculate the Green's function G,^,y of 

CFT resuhs for Gf^fi' and F^^'(e) in the e/o basis; from Eq. (||), and produces the following results for U^^'{e) 

[extracted from Gfjfj' using Eqs. (H3) and (Eq)]: 



these, the above Grj-ri' and Urjrj'{e), which were defined in 
the _L/i?-basis, can be obtained by 



^W' — ^^rjrj^ilV V'v' ' 



(45) 



Ufjfi'i^) has the form 



u' 



u" 



(46) 



In Appendix ^ the above formalism is illustrated by where U " and U ° can be interpreted as the magni- 
a simple example, namely potential scattering of two tudes of the scattering amplitudes in the even and odd 



species of fermions (i.e. r\— 1, 2). 



E. Result of CFT calculation for G^j^j' and Ufjfji{e) 



channels, respectively. Since the odd channels decouple, 
U ° —I. For the even channels, U " has the following 
scaling form: 

C7<^' (e, T) = Ari/2f(£/T)e'*= . (47) 



The CFT calculation of Gfjfj' and Ufjfj'{e) in the e/o Here e'"'^" is a trivial phase shift ^^ that can occur in 
basis, which follows closely the work of AL, is outlined the Kondo channel if particle-hole symmetry is broken 



'^''We shall assume that the phase shift (j}e is energy- 
independent. In general, it can have an energy-dependence. 
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(seetll*'^'^**""^^), and A is a non-universal constant (called 
A7 in paper III). T{x) is a universal scaling function, 
whose explicit form was calculated by ALIHI: 



f(.x) 



3 

2^2 



;27r)i/22sin(7r/2) 



du 



(48) 



^(-.x)/(2.)^^l/2(^_^)l/2j^(^) 



T^3/2) 



-1/2 



(l-u) 



-3/2 



F{u) = F(3/2, 3/2, 1; u) is a hypergeometric function. 
The Jdu integral can be done numerically for any value 
of X, thus giving us an explicit expression for the scal- 
ing function T{x). The real and imaginary parts of 
r(a::) = Te{x) + iTo{x) have the properties 



^e/o{x) =±re/o(-x) , 



and 



Teix) < . (49) 



V. CALCULATION OF THE CURRENT AND 
SCALING FUNCTION 

We now have all the ingredients for step (S6) of sec- 
tion |l^, the actual calculation of the current, from which 
we shall extract the desired scaling function r(a:;). 



A. Calculation of the current 

Using Eq. ( |45| ) to express the current / of Eq. (^) in 

terms of the e/o scattering amplitude Ufjfj'{e) of Eq. (46), 
we find: 



^=^EA^'M^.'(£') + a%,]fis',rj'), (50) 
P^.(£') = X1 (^^^^^) (tT^[n'<'&N^ . (51) 



Note that in spite of the fact that the current operator is 
diagonal in 77 [see Eq. (j^)], -P,, turns out to have o/e cross 

terms, U U . This is a direct consequence of the L-R 
scattering matrix v^^i introduced in Eq. (|^): it necessi- 
tated the L/R-to-e/o basis transformation N^^, and this 
produced a current operator that is off-diagonal in the 
e/o basis. The presence of o/e cross terms in P,y{e') is 
extremely important, since U , describing Kondo scat- 
tering in the even channel, has a T^ " contribution, but 
U , describing no scattering at all in the odd channel, 
does not. Thus we see that our model contains a T^'^ 
contribution to the current, as observed in experiment 
(compare property (Cu.6) in paper I). 

On the other hand, had we attempted to use a model 
without L-R scattering, i.e. with Va-a-' = ^W' (such as the 
model studied by Schiller and HershfieldO), no L/R-to- 
e/o basis transformation would have been needed; then 
P^'(e') would be proportional to U^U, i.e. to (T^/^)^, 
not T^/^. Thus, the inclusion of L-R scattering into the 
model is absolutely essential to obtain the T^/^ depen- 
dence. 

In the above presentation, we glossed over one impor- 
tant subtlety: the scattering matrix Ur/ri' must be unitary 
[see Eq. (|37|)], but the form given in Eq. ([47| ) manifestly 
is not (since, e.g. U " — foif-jT = 0). This reflects the 
so-called "unitarity paradox" EJ, according to which the 
scattering matrix for free fermions off a 2-channel Kondo 
impurity into free fermions is not unitary, which seems to 
violate the conservation of probability dur.inK a scattering 
process. The resolution of this paradoxESEJo is that the 
"missing probability" is scattered into a sector of Hilbert 
space that cannot be described in terms of linear combi- 
nations of single-particle fermionic excitations (compare 
footnote 9) but has a simple representation when the the- 
ory is bosonized. In paper III we shall discuss this issue 
in more detail, and show how to incorporate the resulting 
complications into the present framework. The upshot is 
that the expression (B2h remains valid. 



Let us now analyze the matrix product of Eq. (p^) in- 
dex by index. All matrices are diagonal in a,i, hence 
the sums y^ in P„' are trivial. Next, consider matrix 
multiplication in the index a. Using Eq. ( |ll| ) for Ng^, 
we find 



Pr 



•,>{e')^a'Rc{u^'"\e')u''\e')) 



(52) 



B. Calculation of scaling function V{v) 



We now have gathered all the ingredients to derive 
the sought-after scaling form for the current and con- 
ductance. Inserting Eq. (M) into Eq. (BG) gives 



/ = 



i^4 fde'^ {Rc [c/'°'(£')] + 1} [foie'-eV/2) - /„(e' + eU/2)] 



(53) 



a(0) 



a(i) 



H — —(pi"' -I- . . ., but this will be very weak (since 
e/ef), and only give rise to subleading corrections in the con- 
ductance, i.e. terms of the form {T^^^ /ef)V^i){V/T). 
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where the factor 4 conies from J^a't' ^^'^ ^^'^ surnV^, , written out explicitly, gives the two terms in the last factor. 
Now, the conductance can be written in the forin|^^| 



G 



dl 



dV 



'4h'{i 



Re 



17'=' (e') 



iRe 



Thus, using Eq. (vm and Mt), G reduces to 



G = 2'- 



l-A7oCos0eTi/2r(7iey/T) 



(55) 



Here we introduced a universal scaling function r(w), 
which is defined as follows in terms of the even part Fg 
of the exactly known function T of Eq. (El): 



(56) 



joThiv) = - fdxf, {x + v/2) l-d^foix)] 



where v = eV/T, x = e'/T, fo{x) = l/(e^ + 1). The pos- 
itive constants 7o, 71 are by definition to be chosen such 
that r(w) obeys the normalization conditions [compare 
Eq. (12) of paper I]: 

r(0) = 1 , r(i;) vs. ti2 has slope = 1 as w^ ^ 00 

(57) 

and a minus sign has been included in the definition ( |5q ) 
of r, since f is negative definite [see Eq. (^)]. 

Thus, we have shown that within the present model, 
the conductance obeys^^ the scaling relation 



G{V,T) = Go + BT^'^T{-iiv) 



(58) 



with the universal scaling function r(w), given by 
Eq. (|5^), known exactly. It is plotted as curve 6 in Fig. 0. 
Note that this function is the same as that found in 
Eq. (19) of paper I (for m — 1 there) by a back-of-the- 
envelope calculation. The reason for this agreement is 
that Fe(a;) also turns out to determine the bulk scatter- 
ing rate r~^(e,r) through the relation 

ATi/2fe(e, T) oc 2 (lmE^(e, T) - ImS^(e, 0)) (59) 
= -(r-i(e,T)-r-i(£,0)) , 

where I]^(e,T) i|S-|the retarded bulk electron self-energy 
calculated by ALllir^'''^^-^''^. This a posteriori justifies the 
assumption made in section VFA.2 of paper I, namely 



u''\^e')\+l]{-d,,fo){e'~eV/2) 



(54) 



that the nanoconstriction conductance will be governed 
byT-i(e,T). 

Note that according to the above calculation and 
Eq. (|58|), the slope of the scaling curve [G{V,T) — 
G{Q,T)]/ BT^/"^ seems to be universal, whereas in exper- 
iment it is not [see Fig. 11(a) of paper I]. The reason 
is that in our calculation we assumed that the impu- 
rity sits exactly at the center of the nanoconstriction, 
where the non-equilibrium between L- and i?-movers is 
strongest, and hence feels the full effect of the applied 
voltage. However, as was explained in Section III of pa- 
per I, an impurity not sitting exactly at the center of the 
constriction experiences an effective voltage aiV , where 
the geometrical constant Oi (of order unity) depends on 
the position of the j-th impurity. When summing over all 
contributing impurities, one thus finds expression (21) of 
paper I, which is simply a sum of terms of the form (pq), 
evaluated at slightly different voltages, corresponding to 
different impurity positions in the nanoconstriction. Our 
lack of knowledge about the a^'s forces us to introduce 
another non-universal scaling factor A^ and use the scal- 
ing form 



G(V,T) = Go + BT^'^T{Av) , 



(60) 



when comparing theory with experiment below (see 
Eq. (13) of paper I). When checking in the next section 
whether epxerimental (or numerical) data for G{V, T) 
obeys this relation, we shall plot it in maximally nor- 
malized form (see Section VI B 3 of paper I), i.e. we shall 
plot 



G{V,T) - G{Q,T) 



{Av) 



1/2 



(61) 



with A determined by the requirement that the asymp- 
totic slope of the resulting function be equal to 1 [com- 
pare Eq. (|57|)]. According to Eq. (|6(]|), curves with differ- 
ent T should all collapse onto each other when plotted in 
this way, and the resulting curve should be identical to 
the universal curve F(i;) — 1 vs. v^''^ . 



^^ To see this, use {d^i fo){£' 4- eV/2) — (9j:'/o)(— e' — eV'/2) and then change integration variables, e' 
term of Eq. (m) . Also recall the sign conventions of footnote 4. 

^^Note that consistency with the sign of the experimen- 
tal zero-bias anomaly requires that B = —2e^//iA7o cosine 
mupti be > 0, i.e. Acos(^e < 0. This is in agreement with 
ALB^' ^^°^ who concluded (for the case (/>e = 0) that A < 
in the regime where the Kondo coupling constant is below its 
critical value, Aj^ < A^, i-e. if one flows towards A^, from the 
weak-coupling regime. 



— e' in the second 
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C. Deviations from Scaling 

It should be emphasized that the scahng relation found 
above is only expected to hold for T/T^ <C 1, because it 
is based on keeping just the leading term in an expan- 
sion of Grj-q' in powers of T/T^- If T/T^ is not ^ 1, 
subleading correction terms proportional to (T/T^-)"", 
that have been neglected in the calculation of G^^' , will 
become important. They will give contributions of the 
form (r/T;f)""r„(u), which will cause deviations from 
scaling. In principle, it would be possible to calculate 
the functions r„(w) within our CFT approach, to obtain 



G(y,T)-G'(o,r) 



^A„ (T/T,.)""r„(«) 



(62) 



n=0 



Here, all a„ are positive and thought to be arranged in 
increasing order. These terms arise from all the possible 
irrelevant operators, generated upon the renormalization 
group flow from the weak coupling to the strong coupling 
fixed point. 

If one takes as starting point for the calculation some 
realistic general anisotropic Kondo model which depends 
on a number of different couplings (such as that emerg- 
ing from Zawadowski and coworkers' RG analysis of the 
original TLS-electron interaction, see Appendix O , the 
amplitudes A„ of the subleading terms will depend on 
these couplings and hence be non-universal. For the pur- 
poses of comparing theory with experiment, these con- 
stants would have to be treated as fitting parameters, 
leading to more freedom than one would want for a mean- 
ingful comparision of theory and experiment. Therefore 
we have focussed in this paper only on the n = term, 
with ao = 0, Ao = 1 (by choice of normalization) and 

To{v)^nv). 

On the other hand, for the sake of comparing our CFT 
calculation with other theoretical methods such as the 
non-crossing approximation discussed below, it is use- 
ful to consider the "idealized" situation that the start- 



ing Hamiltonian is the isotropic 2CK model of Eq. (12). 
Then there really is only one coupling constant, vk^ 
which determines T^, and all the amplitudes A„ are ex- 
pected to be universal. This implies that the (maximally 
normalized) differential conductance of (p3) would be a 
universal function of two scaling variables. 



G{V,T)~G{Q,T) 
STi/2 



= r(T/r,.,y/r), 



(63) 



which could in principle be meaningfully compared to 
other calculations. 



VI. FINAL RESULT FOR SCALING CURVE 

In this section we compare the CFT prediction ( |5q ) 
for the universal scaling curve T(v) to the experimen- 
tal scaling curve of Fig. 11(b) of paper I. We also com- 
pare it to the results of Hettler, Kroha and Hershfield 



(HKH)ll3't3, who computed the non-linear conductance 
within the non-crossing-approximation (NCA). 



A. A Few Words on the NCA Method 



In order to understand what HKH did, a few introduc- 
tory remarks about the NCA method and a summary of 
HKH's results are in order here. Some more details may 
be found in Appendix O. 

HKH adopt an SU{2), 2-channel Anderson Hamilto- 
nian in the slave-particle representation [see Eq. (prl]|)] 
that in the infinite-C/ limit that can be mapped by 
a Schrieffer- Wolff transformation onto the NTKM of 
Eq. (H) and (^). The two models are therefore in 
the same universality class and describe the same low- 
energy plaisics. HKH treat their model with the NCA 
techniquetj, which they generalize to 1/ ^^ using 
Keldysh techniques. 

The NCA approach is a self-consistent summation of 
an infinite set of selected diagrams, and hence was con- 
ventionally viewed as an uncontrolled approximation in 
the sense that in general there is no small parame- 
ter. However, when applied to the TV ^ oo, /c — > oo 
limit, with k /N fixed, of the SU{N), fc-channel Anderson 
model (Eq. (HI) with k channels of electrons i — 1, . . . ,k, 
each with N possible pseudo-spin values a — 1, . . . ,N) 
[which maps under a Schrieffer- Wolff transformation onto 
the [/(I) X SU{m,x SUik)f Kondo model], the NCA was 
recently shownE£l to yield exact results for some quantities 
(as discussed in more detail in Appendix O). In partic- 
ular, pJpr all N and k (with fc > 2) the NCA approach 
givea23 the same leading critical exponents as conformal 
field theory (i.e. the exact ones) for all physical properties 
involving the 4-point slave particle correlation functions, 
including the impurity spectral function Ad{Lu) which 
in HKH's calculation governs the point-contact conduc- 
tance. (For example, for general N, k both CFT and the 
NCA yield a — -j^^j^ for the exponent in the scaling re- 
lation (|^), which reduces to ^ for fc = A^ = 2.) The 
NCA method can be therefore be regarded as a useful 
interpolation between the regime T/T^ <C 1, where it 
gives the correct exact critical exponents, and the regime 
T/Tk ^ 1, where any perturbative scheme works. More- 
over, when combined with the Keldysh technique, it deals 
with the non-equilibrium aspects of the problem in a 
more direct way than our CFT approach (it can be re- 
garded as a self-consistent determination of the scattering 
amplitudes), and is able to go beyond the weakly non- 
equilibrium regime (V ^ 2\). 

Therefore, it is certainly meaningful to compare the 
NCA results of HKH to ours. CFT serves as a check on 
how well the NCA does at F = and very low tempera- 
tures, where CFT is exact and NCA only an uncontrolled 
approximation. As we shall see, this check confirms the 
reliability of the NCA method in the regime of very low 
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energies. Conversely, the NCA can then be used as a 
check on our use of CFT for y 7^ situations, for which 
the NCA presumably is the more suitable method. 

Nevertheless, a word of caution is warranted when 
comparing NCA and CFT results: they cannot be ex- 
pected to agree exactly, even deep in the scaling regime 
T -^Tk, since the Anderson model on which the NCA is 
based breaks particle-hole symmetry, whereas the 2CK 
model and its CFT solution does not. Moreover, the 
NCA's deviations from scaling occur when T/Tj^ is no 
longer <^ 1, but this is a regime in which the quality of 
the NCA can not be checked using CFT, and where its re- 
sults for N = k = 2 can well differ somewhat from those 
for the lar ge- A ^, k limit where it can be trusted better. In 
Appendix E5, we suggest a way around these problems, 
based on a new approximation proposed by ParcoUet, 
Georges, Kotliar and Sengupta (PGKS)EZl, which is rem- 
iniscent of the NCA but preserves particle-hole symme- 
try, and moreover becomes rigorously exact in the limit 
A^ ^ cxD, fc ^ 00 at fixed k/N. 

On the other hand, when comparing theory to exper- 
iment, the above caveats are probably less important 
than uncertainties of another kind: as mentioned after 
Eq. (p^), the fact that deviation- from-scaling contribu- 
tions are non-universal means that they will always de- 
pend to some extent on unknown quantities such as the 
initial values of irrelevant couplings. In this light, the 
NCA probably should be viewed simply as a reasonable 
(if non-rigorous) interpolation between to well-checked 
limits, as mentioned above. 



B. NCA Results of Hettler, Kroha and Hershfield 

HKH calculated the conductance G{V, T) for a series 
of temperatures, measured in units of T^, ranging from 
T/Tk = 0.003 to 0.5. Fig. |(a) shows their resuhs for 
G{V,T), plotted according to Eq. (|l|) with A = 1 (i.e. 
without any adjustable parameters). The experimental 
data for one of the Cu samples of Ralph and Buhrman 
(called sample #1 in paper I), for which T^ — 8K, are 
shown for comparison in Fig. 0(b). 

The lowest T/T^ values in Fig. ^(a) show good scaling, 
in accord with the CFT prediction. However, for larger 
T-values, marked deviations from scaling occur, just as 
seen in the experimental curves of Fig. ^(b). It is one 
of the strengths of the NCA method that these devia- 
tions from scaling are automatically obtained, without 
the need for making a systematic expansion in powers 
of T/Tk and V/Tk, as would be necessary in the CFT 
approach. 

The striking qualitative similarity between the two sets 
of curves in Fig. |4| can be made quantitative by using T^ 
as a fitting parameter: the choice of T^ determines which 
curves in Fig. ^(a) and (b) are to be associated with each 
other. Choosing T^^ SK for sample 1, HKH were able 
to get "quite good"lla simultaneous agreement between 



a significant number of the individual experimental data 
curves and their NCA curves of cocEesponding temper- 
ature. This is illustrated in Fig. 013 for 3 curves from 
sample ^ 1. In other words, hy using a single fitting pa- 
rameter, Tfc, HKH obtained good quantitative agreement 
between the NCA and experimental conductance curves 
for a whole set of curves. 



C. Comparison of CFT and NCA Results with 
Experimental Scaling Curve 

Let us denote the result of plotting a given NCA nu- 
merical G{V, T) curve in the maximally normalized form 
of Eq. (|l|) by r{v,T) - 1. Fig. |(a) shows that for suf- 
ficiently small T, the T{v,T) — 1 curves for different T 
all overlap, i.e. the NCA results show good scaling as 
T — > 0, in agreement with the CFT prediction. The 
r{v,T) curve with the smallest T calculated by HKH, 
namely T/T^ = 0.003, is the most likely to agree with 
the CFT result for Viv)] the reason is that for this curve 
the T /Tk deviations from perfect scaling, which are ne- 
glected in the CFT calculation [i.e. T{v) = r(u,0)], are 
smallest. 

In Fig. H we show the three experimental scaling curves 
of Fig. 11 of paper I (curves 1-3), the CFT prediction for 
r(-y) — 1 from Eq. ( pq ) (curve 4), and the NCA result for 
T{v, T) - 1, for T/T^ = 0.003 (curve 5) and T/T,, = 0.08 
(curve 6). All these curves have been rescaled into the 



"maximally normalized form" of Eq. ( |57[ ). We see that 
there is rather good agreement between the CFT curve 
and the T/T^ = 0.003 NCA result. The experimental 
scaling curves agree with neither of these, but agree re- 
markably well with the T/T^ = 0.08 NCA curve. 

To make these statements quantitative, we compare 
the values for the universal constant Fi, defined as fol- 
lows from the asymptotic large-u expansion of T(v) — 1 
[compare Eq. (26) of paper I]: 



r(w)- 1 = 1-1/2 + ri + o(w-i/2). 



(64) 



Fi is the y-intercept of the asymptotic slope of the curve 
r(w) — 1 vs. w^/^, extrapolated back to w = 0. It mea- 
sures "how soon the scaling curve bends up" towards 
linear behavior, and is the single parameter that most 
strongly characterizes the scaling function (which is oth- 
erwise rather featureless). We find the following values 
for Fi: 



r?''^ = -1.14 ±0.10, 

Y^'^^(T/T^ = 0.003) = -1.12 ± 0.10 , 

Ff ^■*(T/rK = 0.08) = -0.74 ± 0.10 , 



Ff^'^ = -0.75 ±0.16. 



(65a) 
(65b) 
(65c) 
(65d) 

Hence, the CFT and NCA calculations for T/T^ = 0.003 
agree rather well, which inspires confidence in the general 
reliability of the NCA method at very low energies. 

The agreement between the experimental curves and 
the T/Tk = 0.08 NCA curve could actually have been 
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anticipated, for the following reason: HKH determined 
their (only) fitting parameter T^ by choosing the value 
(namely T^ = 8 K) that produces the best agreement be- 
tween the few lowest-T curves in their set of calculated 
G(V, T) curves and the corresponding experimental ones. 
Thus, the very lowest T-curve in the experiment (with 
T = 0.6K) is well-reproduced by the corresponding NCA 
curve (with T/T^ — 0.08) because T^ was specifically 
chosen to produce this agreement. ^^ 

It is somewhat surprising, though, that the difference 
between the T/Tj, = 0.003 and T/Tj, = 0.08 NCA curves 
is so large. Perfect scaling would require all the various 
r(w,T) curves for different T to overlap, and the fact 
that they do not shows that the deviations from perfect 
scaling which are expected to develop as T/T^ grows are 
already significant at values as small as T/T^ = 0.08. 

Thus, as first pointed out by HKH, the NCA results 
imply that the T/Tk corrections to the universal scal- 
ing cur ve th at were neglected in the CFT calculation (see 
section VC) are in fact not negligible in the present ex- 
periment: T is still large enough that they matter, and 
the experimental scaling curve is not the truly universal 
one. This conclusion explains why the CFT and exper- 
imental scaling curves don't agree; it also suggests that 
if the experiments were repeated at lower temperatures, 
better agreement might be achievable. 

One might ask whether our conclusion that deviations 
from scaling are important are not in conflict with the 
claims in paper I [property (Cu.6)] that the experimen- 
tal curves show good scaling. The answer is that while 
the experimental curves do scale well, they do not scale 
quite well enough to reproduce "perfect" scaling. Per- 
fect scaling requires that the curves overlap completely 
when plotted in maximally normalized form (as in fig- 
ure 0) , a procedure that involves rescaling the x-axis by 
a constant A to make the slope = I. This procedure 
is clearly very sensitive: even curves that seem to col- 
lapse well onto the same scaling curve when not maxi- 
mally normalized [as those in Fig. ||(b), or Fig. 8(b) of 
paper I] , can show slight differences in slope in the regime 
of largish v when they begin to bend away from the ideal 
scaling curve (note that some uncertainty is involved in 
determining this slope, since the curves are not perfectly 
linear in this regime). When being brought into max- 
imally normalized form, these curves will have their x- 
axes rescaled by different amounts to make all the slopes 
equal to 1 (the exact amount of rescaling needed being 
subject to the same uncertainty as the slope), and can 
by this rescaling be sufficiently deformed that they do 
not collapse onto each other any more. This is vividly il- 



lustrated by the observation that the T/T^ — 0.003 and 
0.08 NCA curves, that in fact seem to overlap rather well 
in the non-maximally normalized form of Fig. H(a), differ 
so markedly in the maximally normalized form of Fig. y. 

In short, maximal normalization is very efficient in re- 
vealing small deviations from perfect scaling, which is 
why the experimental data, which scales well when not 
maximally normalized, does not scale so well under max- 
imal normalization. 

One might be tempted to compare the CFT curve with 
experiment in non-maximally normalized form, where de- 
viations from scaling do not reveal themselves so glar- 
ingly. However, this would not be meaningful, because 
the slope of the CFT scaling T{v) curve is universal, 
whereas those of the experimental scaling curves are not 
(see Fig. 11(a) and the last paragraph of section |6|). 
The only meaningful comparison between CFT and ex- 
periment is in a form in which the non-universality of 
the experimental slopes has been rescaled away, i.e. the 
maximally normalized form. 

^From a theorist's point of view, the conclusion that 
the experimental scaling curve is not the universal one 
and that non-universal T/T^ corrections do play a role is 
somewhat disappointing, since for a system about whose 
microscopic nature so little is known, the quantities that 
allow the most compelling comparison between theory 
and experiment are universal quantities, which are inde- 
pendent of the unknown details. However, disappointing 
or not, this is the message of Fig. 13. 

Nevertheless, the good agreement between the CFT 
and NCA scaling curves, which confirms the NCA 
method to be reasonably reliable [without fully elimi- 
nating the caveats mentioned in the previous section, 
though] , combined with the good quantitative agreement 
between the NCA and the experimental conductance 
curves when T^ is used as fitting parameter, allows the 
main conclusion of this paper: 

The 2-channel Kondo model is in quantitative agreement 
with the experimental scaling GiV, T) data. 



VII. SUMMARY AND CONCLUSIONS 

The calculation of this paper was inspired by experi- 
ments of Ralph and Buhrman on ZBAs in quenched Cu 
nanoconstrictions (reviewed in paper I) . These are quali- 
tatively in accord with the assumption that the anomalies 
are caused by two-level systems in the constriction that 
interact with electrons according to Zawadowski's non- 
magnetic Kondo model, which is believed to renormalize. 



'^'^The NCA calculations achieved more, though, than merely 
fitting one curve with one parameter, because they succeeded 
in repr oducin g quite well a whole set of curves (see end of 
section VIA). 
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at sufficiently low temperatures, to the 2CK model. 

To obtain a quantitative test of this interpretation of 
the experiment, we performed a calculation of the non- 
linear conductance G(V, T) of a nanoconstriction con- 
taining 2-channel Kondo impurities, in the weakly non- 
equilibrium regime (weakly non-equilibrium regime) of 
y, T <C Tk, and extracted from it a certain universal scal- 
ing function T{v), which we compare with experimental 
scaling function. 

To model the experimental situation, we introduced a 
generalization of the bulk 2CK model, namely the nacon- 
striction 2-channel Kondo model (NTKM), which keeps 
track of which lead (left or right) an electron comes from 
and is scattered into. 

The main conceptual challenge in the calculation of 
G(V, T) was how to deal with the non-equilibrium as- 
pects of the problem. On the one hand, standard per- 
turbative Keldysh approaches do not work for T <C T^, 
where perturbation theory breaks down for the Kondo 
problem. On the other hand, Affleck and Ludwig's con- 
formal field theory solution (CFT) of the 2CK problem 
was worked out only for an equilibrium electron system. 

Therefore we proposed a conceptually new strategy 
(outlined in section III, the heart of this paper) which 
combines ideas from CFT with the Hershfield's Y- 
operator formalution of non-equilibrium problems: Her- 
shfield showed that the calculation of non-equilbirium ex- 
pectation values becomes simple when they are expressed 
in terms of the scattering states of the problem. We ex- 
pressed these in terms of certain scattering amplitudes 
Un^Y , which we extracted from an equilibrium two-point 
function G,,,,' — —{Tipriipli') that is exactly known from 
CFT. (This procedure only gives their V = values, but 
we proposed that in the non-Fermi-liquid regime the cor- 
rections of order V/Tj^ are negligible.) Once the Urj^' 
were known, the calculation of the current was straight- 
forward. 

In the present paper, we implemented all parts of this 
strategy, except that which requires a detailed knowledge 
of CFT, namely the calculation of Gj^^' . This is discussed 
in detail in paper III. 

Our result for the scaling curve T(v) does not agree 
with the experimentally measured scaling function, be- 
cause terms of order T/T^ that are neglected in our cal- 
culation are apparently not sufficiently small in the ex- 
periment; however, when our results are combined with, 
the numerical results of Hettler, Kroha and IIershfieldE£l 
(which implicitly do include the neglected terms), quan- 
titative agreement of the 2CK calculations with the ex- 
perimental scaling results is achieved (see section VI C| ). 

Thus we conclude that the NTKM is in quantitative 
agreement with the experimental scaling data. This lends 
further support to the 2CK interpretation of RB's ex- 
periments, and the associated conclusion that they have 
indeed observed non-Fermi-liquid behavior. 

However, the theoretical justification for assuming that 
the non-magnetic Kondo model will under renormaliza- 



tion flow into sufficiently close proximity of the non- 
Fermi-liquid fixed point of the 2CK model has recently 
been called into quesiiniL There are presently unresolved 
theoretical concernsE^lEy, summarized in Appendix y, 
whether a realistic TLS-electron system will ever flow 
sufficiently close to this fixed point to exhibit the associ- 
ated non- Fermi-liquid behavior, because of the inevitable 
presence of various relevant perturbations that can pre- 
vent the flow towards this fixed point. 

Therefore, not all questions regarding the Ralph- 
Buhrman experiments have been resolved to everyone's 
satisfaction. In our opinion, the outstanding question 
that remains is: why does the 2CK interpretation of this 
experiment seem to work so well despite the concerns 
about the theoretical justification for assuming proxim- 
ity to the 2CK model's non-Fermi-liquid fixed point? In 
view of the fact that at present no alternative explana- 
tion for the experiment is known that is in agreement 
with all experimental facts, we believe that the question 
of the flow towards and stability of the non- Fermi-liquid 
fixed point of the non-magnetic Kondo problem is worthy 
of further theoretical investigation. 
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APPENDIX A: SEMICLASSICAL DESCRIPTION 
OF NON-EQUILIBRIUM TRANSPORT 

In order to motivate th e for m of the free density matrix 
Po introduced in section 11 B, we recaU in this appendix 
some standard results from the semi-classical theory of 
non-equilibrium transport through a ballistic nanocon- 
striction. Usually, this is described using a semi-classical 
Boltzmann formalism to calculate the semi-classical elec- 
tron distribution function f^{r) and the electrostaJic_po- 



tential energy e(f>{r). This was first worked out iiiE3'E2l; a 
very carefuL4|reatment may be found iiiE3, jwhich is well- 
reviewed inEj. A more up-to-date review iaij. 

In the semi-classical strategy, one first calculates 
/i (r) and 60'-°-' (r), the distribution function and elec- 
trostatic potential in the absence of any electron scat- 
tering mechanism, and thereupon uses these functions to 
calculate the backscattering current due to electrons that 
are backscattered while attempting to traverse the hole. 
The results for /i ' {r) and e(f>''^\r) are standard and 
shown in Figs. and 0.^* Fig. pi is a position- momentum 
space hybrid, showing f-{r) dX T — 0, with its fc-space 
origin drawn at the position r to which it corresponds. 
One can understand Fig. ^(a) almost without calculation, 
simply by realizing that in the absence of collisions, elec- 
trons will maintain a constant total energy E^. Thus, 
an electron that is injected from z = ±cxd in the R/ L 
lead with total energy £^^(z = ±oo) ~ e^± eV 12 and 
traverses the hole, will experience a change in its poten- 
tial energy from e<^(±oo) = ±6^/2 to e0(=Foo) = '^eV jl 
and hence accelerate or decelerate in such a way that 
E^{r) = Eg + e(j){f) remains constant. 

The key feature of Fig. is that the distribution of 
occupied electron states in momentum space, at any 
point r in the vicinity of a ballistic constriction, is highly 
anisotropic and consists of two sectors, to be denoted 
by L and R. The L/ R sector contains the momenta of 



all electrons that are incident as L/ R-movers, i.e. origi- 
nate from the ±y/2 or R/ L side of the device, and have 
reached r along ballistic straight-line paths, including 
paths that traverse the hole (the bending of paths due 
to the electric field is of order eV/ep and hence negligi- 
ble). At a given point r, the momentum states in the L/ R 
sectors are filled up to a maximum energy of (EAr)] 
which, because of energy conservation along trajectories, 
is equal to ^ ± 6^/2, the Fermi energy at z = ±oo from 
where the electrons where injected. Thus, for k in the 
L/R sector, one finds 

= /o [e^ - (// ± eV/2 - e(/.(°) (f)) ] . (Al) 

Fig. H shows that et/)'-*'-' (r) changes smoothly from ~eV/2 
to -\-eV/2 (the change occurs within a few constriction 
radii a from the hole). It is worth emphasizing, though, 
that the electrostatic potential energy e(j){f) plays only 
an indirect role when it comes to calculating low-energy 
(i.e. T /eF,V/£F <C 1) transport properties. The reason 
is simply that the only role of e0(r) is to define the bottom, 
of the conduction band, hence causing acceleration and 
deceleration of electrons to maintain E^{r) — constant. 
Low-energy transport properties, however, are deter- 
mined by what happens at the top of the conduction 
band, in particular by the sharply anisotropic features 
characterizing Fig. and Eq. (Al).^^ 

The above considerations suggest that the essence of 
the non-equilibrium nature of the problem will be cap- 
tured correctly if we adopt the following simplified pic- 
ture: ignore the spatial variation of the electrostatic po- 
tential e(f)(r) altogether, and simply consider two leads 
{R/L) with chemical potentials (measured relative to the 
equilibrium /i) /i^ = a^eV, which inject i/i?-moving 
ballistic electrons into each other (recall that a = (-I-, — ) 
for (i, i?)-movers). The two leads are assumed infinitely 
large and hence "independent and unperturbed" , in the 



'^^Our figures and arguments are given for the case eV > 0. 
We take e = -|e| and hence V = -\V\. With fi ± eV/2 for 
R/L leads, there then is a net flow of electrons from right to 
left, and the current to the right is positive. 

^^This is illustrated, for example, in the calculation of the 
Sharvin formula for the conductance Go of the sr-s:ircular con- 
striction (radius a) in the absence of scatteringEj: 



/o = / dxdy — ^(«p)z/,?(. 

<J hole ^ 



x,y,z = 0) = a e msF/{2nh )\V\ 



(A2) 



It depends on the electrostatic potential only through 

e(j>{x,y,z = 0) = 0, and it is easy to verify that the 

y-dependence arises solely from the L/R anisotropy of 
fpix,y,z = 0). 
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sense that their thermal distribution properties are not 
perturbed when a sniaU number of electrons are trans- 
ferred from one to the other. This simplified picture is 
the basis for the Ansatz (||) for the free density matrix 
Po in section II B. 



APPENDIX B: THE BULK NON-MAGNETIC 
KONDO MODEL 

In this appendix we recall some basic properties of 
Zawadowksi's non-magnetic (or orbital) Kondo model 
for the interaction of a TLS with conduction electrons 
in a bulk metal. Thus, this appendix provides the 
background material that was assumed known when 
we introduced the nanoconstriction two-channel Kondo 
model [in section y. Zawadowski proposed his model 
in RefM— suhsfioujently developed it with his coworkers 
in Refs.olljOEj, and rather recently, together .-witii 
Zarand, introduced some important refinementsE3Ea. 
I and exhaustive reviews may be found 
a, respectively. 



1. Zawadowski's Bulk Bare Model 

Consider a tunneling center (TC) in a bulk smetal, 
i.e. an atom or group of atoms that can hop between 
two different positions inside the metal, modelled by a 
double- well potential [see Fig. |l|, and Fig. 6 of paper I]. 
At low enough temperatures and if the barrier is suf- 
ficiently high, hopping over the barrier through thermal 
activation becomes negligible. However, if the separation 
between the wells is sufficiently small, the atom can still 
move between them by tunneling. 

If the tunneling is slow (hopping ratesLJ r^^ < 
10®s~^), the atom is coupled only to the density fluc- 
tuations of the electron. sea, which can be described by 
a bosonic heat batl£3C3. The tunneling is then mainly 
incoherent, and the only effect of the electron bath is 
to "screen" the tunneling center: an electron screening 
cloud builds up around the center and moves adiabati- 
cally with it, which leads to a reduced tunneling rate due 
to the non-perfect overlap of the two screening clouds 
corresponding to the two positions of the tunneling cen- 
ter. 

In this paper we are interestecL-Qnly in the case where 
the tunneling is fast (at ratea^ lO^s^^ < t~^ < 
lO^^s^^), in which case the tunneling center is usually 
called a two-level system (TLS) [though in this appendix 
and the next we shall continue to call it a tunneling cen- 
ter, because in general m ore than two states can be as- 
sociated with it, see Eq. (|B^)]. Then the energy corre- 
sponding to the tunneling rate, determined by the un- 
certainty principle, is in the range 1 mK to 10 K. (If 
the tunneling is "ultra-fast" (t~^ > lQ^'^s~^)^ the energy 
splitting E2 — El between the lowest two eigenstates due 



to tunneling becomes too large (> IQK) and the interest- 
ing dynamics is frozen out.) Moreover, the TLS-electron 
coupling is assumed strong enough that in addition to 
screening, an electron scattering off the tunneling cen- 
ter can directly induce transitions between the wells: it 
can either induce direct tunneling through the barrier 
{electron-assisted tunneling)^ or excite the atom to an 
excited state in one well, from where it can decay across 
to the other well {electron- assisted hopping over the bar- 
rier) . 

To describe such a system, Zawadowski introduced the 
following model. The Hamiltonian is the sum of three 
terms: 



H = Hrr 



H, 



i?.„* . 



(Bl) 



The first term describes the motion of the tunneling cen- 
ter the double well, in the absence of electrons [see Fig. 1^, 
and Fig. 6 of paper I]: 



H^c = J2 ^'^^^^^ 



(B2) 



This problem is considered to be already solved: the en- 
ergies Ea {El < E2 < . . .) are the exact eigenenergies of 
the exact eigenstates |^o) = ^I|0) of the tunneling center, 
with corresponding wave- functions ipa{R)- The spectrum 
will contain two nearly-degenerate energies Ei and E2, 
split by an amount E2 — Ei, corresponding to even and 
odd linear combinations of the lowest-lying eigenstates 
of each separate well; the remaining energies, collectively 
denoted by E^^ , correspond to more highly excited states 
in the well, with E^^ — E2 typically on the order of the De- 
bye temperature of the metal, i.e. several hundred Kelvin. 
The tunncling-center-electron interaction is described 
by a pseudo-potential V{R — r), which describes the 
change in energy of the tunneling center at position R 
due to the presence of an electron at position r, and is 
assumed to depend only on the relative coordinate r — R: 

H,nt = J2 fd.''r¥d^■^^{r)V{f- R) 

X ldRY,blKiIi)ba'Va'{R) , (B3) 

aa' 

where ^i(af) = (Vol)^^/^ ^-.e'^'^Cop'i . Here c]^-.. creates 
a free electron (hence the subscript o) with momentum p 
{—pp), energy Sp (assumed independent of the direction 
p or rip of p) and Pauli spin i =t, i (we use the index 
i because this will turn out to be the channel index). 
Terms in Eq. (|B3| ) with a ^ a' correspond to transitions 
between eigenstates of the tunneling center induced by 
the scattering of an electron. 

Now, let {Fa{p)} be any complete set of orthogonal 
functions of p (e.g. Fa{p) = V4TrYim{p), but in princi- 
ple any set of orthogonal angular functions can be used) , 
labelled by a discrete index a and satisfying 
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/ 






(B4) 



Then the electrons' continuous direction index p can be 
traded for the discrete index a by making a unitary trans- 
formation {No is the density of states per spin at Sp): 



Co... = N^' J ^ F:{p)c 
Cop. = N-^/^J2F^{p)co, 



(B5) 



Written in this form, the interaction has the form of a 
generahzed, anisotropic Kondo interaction: a and a can 
be regarded as impurity- and electron pseudo-spin indices 
(since a takes on infinitely many values, the electrons 
have an infinitely large pseudo-spin), and the interac- 
tion describes electron-induced "spin-flip" transitions of 
the impurity. Note, however, that because the nature of 
the interaction is non-magnetic (to which fact the model 
owes its name), the interaction is diagonal in the Pauli 
spin index i = (t,i)- Thus we have two identical chan- 
nels of conduction electrons, the {c^^jj|}- and the {c^j^i } 
operators, and accordingly i is called the channel index. 



The new set of operators {coeai} are labelled by the con- 
tinuous energy index e (= £p) and the discrete index a, 
to be called the conduction electron pseudo-spin index, 
for reasons that will become clear below. 

In the new basis, the electrons' kinetic energy and in- 
teraction with the tunneling center can be written in the 
following form: 



H, 



Hi. 



J -D 



'I aa' aa' 



^ c 

D I'D 

df I df' v°-°-',r^ c , ,bH , 
-D J-D 



(B6) 

(B7) 



For simplicity, the standard assumptions were made that 
electron energies lie within a band of width 2D, symmet- 
ric about Ef, with constant density of states No per a,i 
species, and that the energy dependence of the coupling 
constants f J|J^, can be neglected. (These assumptions are 
justified by the fact that the Kondo physics to be stud- 
ied below is dominated by excitations close to the Fermi 
surface.) The w°°, are volume-independent, dimension- 
less constants (typically of order 0.1 or smaller), whose 
actual values are determined by the potential V{r — R) 
and the tunneling center eigenstates ipa{R)- 



sas' a' ^aa' 



E 

b/3 



D 
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2. The Renormalized Bulk Model 



The formal similarity of the interaction of Eq. (B7) 
with the usual (magnetic) Kondo interaction implies that 
here too perturbation theory will fail at temperatures be- 
low a characteristic Kondo temperature T^, leading to 
complicated many-body physics as T -^ and a strongly 
correlated ground state. Perturbation theory fails for 
T < Tk because the effective (T-dependent) coupling 
constants u°°, grow as T decreases, and eventually be- 
come too large (see Fig. H in Appendix O). The way in 
which this happens was studied in great detail by Za- 
wadowski and co-workers. Using Anderson's poor man's 
scaling technique to analyse the renormalization group 
evolution of the bare model, they concluded that the 
renormalizecLaiiodel to which it flows as the temperature 
is lowere d isE I the isotropic two- channel Kondo model 
(see Eq. ( |B10| ) below) . Below we briefly give the starting 
point and final result of their poor man's scaling analysis. 
A summary of the intermediate steps and main assump- 
tions made along the way can be found in Appendix y. 

The interaction vertex, calculated to second order in 
perturbation theory, is given by the following expression: 



"^ ''"'£' + K'-(e+i?fc) 






Ip{e) 



/3a' ^aP 77 



e' + Ea' - {-£+e'+e+Eb) 



^ln[max{K',S6,r,£,£'}/i?] 

b/3 



afc ha' _ „,ba' ab 



(B8) 
(B9) 



To obtain the second line, only the logarithmic terms 
were retained. 



Note 

,,ab ^,ba' 



the 



^ap'"l3a' ''^al3 ^ 0a' 



occurence 

ah 



of the "commutator" 
the fact that this is in general 

non-zero, due to the non-trivial angular dependence of 
the coupling constants, is crucial for the presence of log- 
arithmic corrections (and is the reason why this model 
is sometimes called a non- commutative model) . 

Now Anderson's poor man's scaling RGE^ is imple- 
mented: one changes the bandwidth from D to a 
slightly smaller D' , and compensates this change by in- 
troducing new coupling constants that depend on a; = 



InD/D', namely v — v{x), with the x-dependence 
chosen such that r°°, remains invariant. The proce- 
dure is repeated successively until D' reaches E^ax — 
niax{Ea' , Ei,,T,e,e'}, at which point the RG flow is cut 
off, and the resulting renormalized model, with coupling 
constants v(hi(D / Emax)) , has to be analyzed anew. 



The upshot of a lengthy analysis (summarized in Ap- 
pendix |c|) is the following: All but the lowest two 
of the excited states of the tunneling center decouple 
from the interaction, which eventually involves an im- 
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purity with effectively only two states, ^"^ a = 1,2 (i.e. 
a TLS with an effective pseudo-spin Simp — 1/2), with 
a renornialized splitting E2 ^ Ei = A and Hamiltonian 
Htc ^ Ea,a'=i,2^1(5^<a')&a'. Likcwisc, for the Con- 
duction electrons all but two of the pseudo-spin degrees 
of freedom, which we label by a ~ 1,2, decouple from 
the interaction. These two 



channels, c! 



survivnig ^..a....i^.^, o^^^j 
and CQ^2i) ^^^ ill general two complicated linear combi- 
nations of the initial c^^^^^s. They represent those two 
angular degrees of freedom that initially were coupled 
most strongly to the impurity and for which the couplings 
hence grow faster under the renormalization group than 
those of all other channels (which hence effectively de- 
couple). Furthermore, the resulting effective interaction 
is spin-isotropic (spin-aaisotropy can be shown to be an 
irrelevant perturbationE3) , so that the e ffecti ve interac- 
tion can be written in the form [see Eq. (C15)]: 



•' •' aQ' = l,2aa' = l,2i=T,i 

^'^K \^oeai\^aa'Coe'a'i\- ip\\Saa'ba') ■ (BIO) 

Here v^ is the magnitude of the effective tunneling- 
center-electron coupling (and estimatedeS'"^"''''^^ to be of 
order v^ ~ 0.1 — 0.2). Thus, the effective Hamiltonian^^ 
has exactly the form of the isotropic, magnetic 2CK prob- 
lem, with impurity pseudo-spin Simp — 1/2 (a = 1,2), 
electron pseudo-spin s^i = 1/2 (a — 1,2), and the Pauli 
spin i =1, 1 as channel index. 

When the temperature is lowered even further, then, 
provided that A = 0, this model flows towards a non- 
trivial, non- Fermi-liquid fixed point at T — Aii-af which 
the system shows non-Fermi-liquid behaviortilO. How- 
ever, A is a relevant perturbation (with scaling dimension 
— ^, see section VI C of I). This means that if A 7^ 0, the 
flow towards the non-Fermi-liquid fixed point will be cut 
off when T becomes smaller than A'^/Tj^, after which 
the flow will be towards a different, Fermi- liquid fixed 
point that corresponds to potential scattering off a static 
impurity. In subsequent sections we shall always adopt 
assumption (A2) of paper I (for reasons explained in Sec- 
tion VI C of paper I) namely that A is sufficiently small 
relative to T (A ^ ^/TT^) that the physics is governed 
by the non-Fermi-liquid fixed point, and that the depar- 
ture of the flow from the latter towards the Fermi-liquid 



fixed point has not yet started. 

It is temp ting to propose for the effective Hamiltonian 
of Eq. (BIO) the following physical interpretation (which 
is given in this form by Moustakas and Fishero, and can 
be viewed as complimentary to Zawadowski's picture of 
electron-induced tunneling). A charged impurity in a 
metal will be screened by a screening cloud of electrons, 
which can be thought of as part of the "dressed" impurity. 
If the impurity is a two-state system, it will drag along 
its tightly bound screening cloud as it tunnels between 
the wells. In doing so, it will redistribute the low-energy 
excitations near the Fermi surface. In particular, it will 
likely interact most strongly with two spherical waves of 
low-energy electrons, "centered" |J3ji the two impurity po- 
sitions in the left and right wellalJ''''^^^^, with which one 
can associate a pseudo-spin index a = L, R. Now, when 
the impurity and its screening cloud tunnels from the left 
to the right well, low-energy electrons around the right 
well will move in the opposite direction to the left well, 
to compensate the movement of electronic charge bound 
up in the screening cloud, and thereby to decrease the 
orthogonality between the pre- and post-hop configura- 
tions. Thus, a flip in the impurity pseudo-spin is al- 
ways accompanied by a flip in electron-^seudo-spin, as in 
Eq. (BIO). In two very recent paperg23, Moustakas and 
Fisher have used this interpretation as a starting point 
for a related but not quite equivalent description of the 
tunneling center-electron systemc3. 



APPENDIX C: POOR MAN'S SCALING 

ANALYSIS OF BULK NON-MAGNETIC KONDO 

MODEL 



In this appendix, we summarize_fpllowing the rergnt 
papers by Zarand and ZawadowskicZlcj and ZarandoLO, 
the poor man's scaling arguments that suggest that the 
bare, equilibrium non-magnetic Kondo model of Eq. (B3) 
renormalizes to the isotropic 2-channel Kondo model 
Eq. ( BIO ). It should be mentioned at the outset, though. 



that the ensuing analysis has a somewhat heuristic char- 
acter, since it employs scaling equations derived in the 
weak-coupling limit, based on perturbation theory in the 
coupling constants. Since such scaling equations cease to 
be strictly valid as soon as one scales into strong-coupling 



■^"The two states are considered here in the energy repre- 
sentation, i.e. their wavefunctions are 931,2 = -js{<Pr i^^pi) in 
terms of the wavefunctions ipr or tpi describing the tunneling 
center locahzed mainly in the r or I wells. 

^'^Of course, the flow toward the isotropic 2CK model only 
happens provided that all relevant perturbations that would 
drive the system away from this flxed point are sufficiently 
small - this implicit assumption will be critically discussed in 
Appendix O. 
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regions of parameter space, by such an analysis the con- 
clusion that the bare model flows towards the 2CK model 
can at best be made plausible, and never be proven con- 
clusively. In fa c t , I this conclusion has recently been called 
into questionll3t3, on the basis-of theoretical considera- 
tions (controversial themselveaHj), that are discussed in 
Appendix O. 



1. Hamiltonian and Initial Parameters 



[this is because the overlap integrals in Eq. ( p^ ) are larger 
for (^*^(/Jr (;■) than for (/3*(/Ji, since the excited state wave- 
function spreads over both wells (see Fig. Q)]. Although 
the amplitudes for such processes are proportional to the 
factor l/£'ex (which is small, since E^^ is large), Zarand 
and Zawadowski showed that such terms also grow un- 
der scaling [see Eq. (p3) below] , and eventually lead to a 
renormalized model which has sufficiently large effective 
tunneling amplitudes to display Kondo physics. 



The starting point is the Hamiltonian introduced 
and motivated in Appendix ^, written in the form of 
Eqs. (^, (^) and (^. Let Ab = E2- Ei be the bare 
energy difference between the two lowest-lying, nearly 
degenerate eigenstates states of the well. The remaining 
energies, Ea, a = 3,4,..., collectively denoted by E,^, 
correspond to more highly excited states in the well. 

We are interested in the regime where A;, <C T ^ 
E^-, ^ D. Hence we take Af, ~ 0, i.e. consider a sym- 
metrical double well with a two-fold degenerate ground 
state. It is then convenient to make a change of basis 
from the exact symmetrical and anti-symmetrical ground 
states 1^1) and [^2) to the right and left states \r) and 

10 =^(|*l) ±1^1/2)). 

Note that, since a non-zero bare tunneling matrix ele- 
ment (Aq) between the wells always leads to a splitting 
El — E2 — Ao, by taking Afc ~ we are also implicitly 
assuming that Ao <C T. This means that direct tun- 
neling events are very unlikely, raising the question of 
whether Kondo-physics will occur at all.^^ However, the 
inclusion of excited states iiL-the model overcomes this 
potential problem as £ollowsoc3: a careful estimate of 
t he coupling constantsE3 in terms of the overlap integrals 
(B3) has shown that 



^.l.ex I 



10" 



-3|„M_„r-,r| 



\^r,ex^ _ , M 



(Cla) 
(Clb) 



The first relation reflects the fact that direct electron- 
assisted tunneling, parameterized by \v^ \, is propor- 
tional to the bare tunneling rate Aq and hence very small. 
However, the second relation shows that the matrix ele- 
ments for electron-assisted transitions to excited states, 
parametrized by |v''''^| and |w''''^^|, are of the same order 
of magnitude as for the usual "screening term" |v''' — w'"''"! 



2. Poor Man's Scaling RG 

The interaction vertex, calculated to second order in 
perturbation theory from the diagrams in Fig. 0(a), is 
given by Eq. (B9): 



eae' a.' aa' 



^ In [max{£;,-, Sfc, T, e, e'}/D] 

bl3 



^ab ^ha' ^.^ci'^.a^ 



V^rV 



aP^pa' ~ "al3^l3a' 



V^rVr 



(C2) 



Now Anderson's poor man's jSJcaling RGE3 is imple- 
mented (very nicely explained \T]^,secUonsS,.2.2y^ electron 
or hole excitations with large energy values do not di- 
rectly participate in real physical processes; their only 
effect occurs through virtual excitions of the low-energy 
states to intermediate high-energy states. Hence such 
processes may be taken into account by introducing 
renormalized coupling parameters, which sum up all the 
virtual processes between a new, slightly smaller cut-off 
D' and the original D. In other words, all virtual pro- 
cesses between the energies D' and D are integrated 
out and their contributions are incorportatcd in new, 
Z?'-dependent coupling constants. This procedure is re- 
peated for smaller and smaller D', until D' becomes on 
the order of maxji^c, T, Spi}. 

Concretely, this is done by writing u°^, = i'aQ'(a;), 
where x = ln{D' /D) and the x-dependence of the cou- 
pling constants is determined by the requirement that the 
interaction vertex be invariant under poor man's scaling, 
i.e. d^r^''^, = 0. By Eq. (^), this leads to the following 
leading-order scaling equation: 



d.lfix) 



Y^6{D' 



E,)[v^^{x),v^\x)] 



(C3) 



■^^This was a serious limitation of Zawadowski's original 
model, which did not include excited states: to give non- 
trivial many-body physic (i.e. a sufficiently large Kondo en- 
ergy Tk), the bare tunneling rate Ao could not be too small; 
yet at the same time, the model only flows to the interest- 
ing non-Fermi liquid fixed point if E\ — E2 <Si T. This would 
have required a rather delicate and perhaps questionable fine- 
tuning of parameters. This problem-iaa been overcome by in- 
cluding excited states in the modelcUca, as explained above. 
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ab 



,,ab 



where we have adopted the matrix notation v^ 

[The significance of6{D' — Ec) is explained in section C 5 



This equation, to be solved with the boundary condition 
w°''(0) — («"'') bare, determines the nature of the RG flow 
away from the weak-coupling limit. 

In the following two sections we outline the results 
obtained by Zawadowski and co-workers concerning the 
nature of the fixed point that the Hamiltonian flows to- 
wards as it scales out of the weak-coupling region. How- 
ever, the arguments that are to follow all have a some- 
what heuristic character: since they are based on scal- 
ing equations that were derived in the weak-coupling 
limit, based on perturbation theory in the coupling con- 
stants, in principle they cease to be strictly valid as soon 
as one scales into strong-coupling regions of parameter 
space. (The only method that gives quantitatively reli- 
able result s far .the cross-over region is Wilson's numerical 
NRGEJT^^MeI.) Many of the results obtained below are 
therefore of mainly qualitative value, and not expected 
to be quantitatively accurate. 



3. Scaling to 2-D Subspace 

Let us for the moment consider the model without any 
excited tunneling center states, i.e. with J2c ~ X^r i (^® 
was done in the first pa persa oM) , postponing the more 

In this case, the coupling 



general case to section C5, 
constants v"'''{x) can be expanded in terms of Pauli ma- 
trices in the 2-dimensional space of the tunneling center 
(in the l-r basis), 



,,ab 



(x) 



A=0 



Vaa'{^)<^ab , 



(C4) 



where A — (0,1,2,3) — {0,x,y,z) and (t';^!^ = Sab- The 
v^ term is called the screening term, and characterizes 
the difference in scattering amplitudes for processes in 
which an electron scatters from the atom in the right or 
left well without inducing a transition to the other well. 
The v^ and v^ terms are called electron- assisted tunnel- 
ing terms, and describe the amplitude for processes in 
which the scattering of an electron induces the tunneling 
center to make a transition to the other well. According 
to Eq. (CI), u"^(0) ~ F(0) < F(0). If one chooses the 
wave-functions of the tunneling center tOr-be real, time- 
reversal invariance requires ^^(0) = (seell3^('')''^'J(^-i^'). 
The problem is now formally analogous to a (very 
anisotropic) magnetic Kondo problem in which a spin- 
i impurity is coupled to conduction electrons with very 
large pseudo-spin (since a takes on a large number of 
valuesl-| However, Vladar and Zawadowski (VZ) have 
shownB('')'^'='=**°"-f"C that (with reahstic choices of the 
initial parameters) the problem always scales to a 2- 
dimensional subspace in the electron's a-index, so that 
the electrons have pseudo-spin Se = h (this happens in- 
depedent of the signs of the initial couplings, see sec- 
tion |C4|). Their argument goes as follows: 



In the notation of Eq. ( C4 ) , the scaling equation ( |C3| ) 
takes the fornMpi"3,e9-(feT 



dx 



= -2j^e^^^t;««^ 

BC 



(C5) 



Since %f{Q) < v^(0) and i^(0) = 0, Eq. (|C|) can 
be linearized in y^ and yy . VZ solved the linearized 
equations in a basis in a-space in which w^n(O) is di- 
agonal Ef a(0) = (5q/3W^(0)], and obtained the following 

S0luti0lMpl"6,eg.3.17. 

Kpi^) = S^pvli^) (C6) 

Kp (^) = <p (0) cosh 2x [vl {Q)-vl (0)] , (C7) 

€p (^) = ^Kp (0) sinh 2x [vl^ (0) - vl (0)] . (C8) 

Barring unforeseen degeneracies in the matrix v_^ , this 
shows that the two elements of y^ which produce the 
largest difference \v^{0) — Va{^)\ will generate the most 
rapid growth in the corresponding couplings v^a{x) and 
^^g(x). In fact, since this growth is exponentially fast, 
any couplings with only slightly smaller |w5(0) — ^^(O)! 
will grow much slower and hence decouple. Thus, we con- 
clude that according to the leading-order scaling equa- 
tions, the system always renormalizes to a 2-D subspace 
in which the electrons have pseudo-spin Se= 2- 

The argument just presented is not quite waterproof, 
though. Firstly, it depends on the assumption of ex- 
treme initial anisotropy in the couplings, and secondly, it 
is based only on the leading-order scaling equations. As 
one scales towards larger couplings, sub- leading terms in 
the scaling equations can conceivably become important. 
Zarand has investigated this issue by including next-to- 
leading-order logarithmic terms [generated by the dia- 
grams iUpEjg. 0(b)] in the scaling equations, which turn 
outtobeEl^'=«P-6): 



d^.yf' = -2i^e 



ABC^B^C 



BC 



~2Nf Y^ {y^TT[{y^f] - v^Triit^v^]) . (C9) 

Note that the number of channels, Nf (equal to 2 for 
the case of interest), shows up here for the first time in 
the next-to-leading order, since each electron loop [see 
Fig. 0(b)] carries a factor Nf. Performing a careful anal- 
ysis of the stability of the various fixed points that oc- 
cur, he concluded that the above-mentioned Se = ^ fixed 
point is the only stable fixed point in of these equations. 
Zarand and Vladar also investigated the effect of the 
other channels, that don't couple as_strongly as the two 
dominant ones, near the fixed pointEll. They produce ir- 
relevant operators that eventually scale to zero (which is 
why these channels decouple), but that can nevertheless 
infiuence the critical behavior near the fixed point. How- 
ever, Zarand and Vladar found that they have the same 
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critical exponent as the leading irrelevant operator in the 
pure 2CK model, which means that these extra operators 
don't change the universal critical behavior, merely some 
of the corresponding amplitudes. 

Since these results are independent of the value of 
Nf and the number of orbital channels considered, and 
Zarand's analysis is exact in the limit Nf — > oo, he ex- 
pects his results to also be valid for Nf = 2. (How- 
ever, no completely rigorous proof exists yet for this ex- 
pectation; in particular, his analysis assumes A;, = 0, 
and-the case Af, ^ is substantially more complicated, 
ggg|io|,(b),section/77 -p^j. anotber, symmetry-based argu- 
ment in favor of S'e = 5, seeEl''*^'=**°"3.3.2(«i) _■) 



4. The fixed point is Pseudo-Spin Isotropic 



Next one shows, foUowin, 



IsecUonlll^ that the Se = ^ 

fixed point is actually isot ropi c in pseudo-spin space. 

The last term in Eq. ( |C9| ) can be eliminated from 
the fixed-point analysis by making a suitable orthogonal 
transformation ir^ -^ TIb^abvP- Therefore, it is suf- 
ficient to c onsid er the first two terms on the right-hand 
side of Eq. (C9). At the fixed point, where dxV^ = 0, we 
have 

Y^ e^^^vj^v^ = ^NfV^ Y^ T,[{v^f] . (CIO) 

BC B=jtA 

Multiplying by w'^ and taking the trace, one obtains the 
three relations 



iNfCx'^ia." + 0.'-') =13 , 

where {A, B, C} = {x, y, z} (cyclically) 



(Cll) 



where we have defined cx.^ = Tr[(w'*)^] and /3 = 
Ti(v_ v^vf^ — vpvPv_ ). This immediately implies one of 
two possibilities: either at least two of the a'^'s are zero, 
which is the trivial (commutative) case without electron- 



assisted tunneling {if 
equal: 



0); or else they are all 



(C12) 



The latter case is the one of present interest. The con- 
clusion that the couplings are all equal (i.e. the effec- 
tive Haiailtonian isotropic) was checked numerically by 
ZarandcJ'^*^'', and is illustrated in Fig. @. 

What is the matrix structure of the 1; s? Introducing 
the notation J^ = tti^w^, Eqs. ( |C10| ) and ( |cT2| ) imply 



2Nfa- 



that the J'^ satisfy the SU{2) Lie algebra. 



(C13) 



which means that they must be a direct sum of irre- 
ducible SU{2) representations: 



J' 






(C14) 



®k=l 



According to the analysis of Zarand mentioned in the 
previous section, only a single subspace 5*6 = 5 in this 
sum corresponds to a stable fixed point (all the others 
correspond to unstable fixed points), in the vicinity of 
which we can therefore write J^^, = \o'aa'- 

After a rotation in a-space to line up the quantization 
axis of the pseudospins of the impu rity and the electrons, 
the effective Hamiltonian to which (P3) renormalizes can 
be written as: 

H,.. ^ fde fde' Y Y Y 

•' '' aa'=l,2aa' = l,2i=],i 

X ""k Wosaih^aa'Coe'a'ij ' {b\\aaa'ba') ■ (C15) 

Here Vk is the magnitude of the effective tunneling 
center-electron coupling (and estimatedcfl'"^''''''^"'^ to be of 
order v^ ~ 0.1 — 0.2). This is the main result of the RG 
analysis: The effective Hamiltonian has exactly the form 
of the isotropic, magnetic 2-channel Kondo problem; the 
two surviving orbital indices a = 1,2 play the role of 
pseudo-spin indices and the Pauli spin indices i =t, i the 
role of channel indices. 

An intuitively appealing interpretation of this model, 
due to M ousta kas and Fisher, is given in Appendix H, 
after Eq. ( |bTo|) . 

We conclude this appendix with a number of miscella- 
neous comments: 

The fact that one always scales towards an isotropic 
effective Hamiltonian is rather remarkable (though in ac- 
cord with the conformal field theory resultS|-t|hat show 
that anisotropy is an irrelevant perturbationEj^'^*''^'^^'): 
the initial extreme anisotropy of the couplings is dynam- 
ically removed, and a SU(2) symmetry emerges that is 
not present in the original problem! 

Note that the initial signs of the anisotropic coupling 
constants did notjjnatter in the above arguments. A more 
careful argumentO''"''^*'""'^'^'^'"^ shows that the flow to- 
ward this fixed point indeed occurs irrespective of the 
initial signs of the coupling constants. 

Relevant perturbations: When the initial splitting 
Afc is non-zero, the r^-nd order RG is consider- 
ably more comphcatedtj''''-'^*'^^*'""'^'^^. The result is 
that Aft gets normalized downward by about two or- 
ders of magnitudetj^'''''^'^'^. However, as emphasized 
jjj52,section3.4.i(c)^ ^/jg splitting Af, is nevertheless a rel- 
evant perturbation: it can be shown to scale downward 
much slower than the bandwidth D' , so that Ab{D')/D' 
grows as D' is lowered. 

By analyzing the stability of the fixed point equations 
against a perturbation that breaks channel symmetry, it 
can likewise be shown that channel anisotropy is a rele- 
vant perturbationEa^'"='=**°"3''-i('=' . 

Kondo temperature: The Kondo temperature is the 
cross-over temperature at which the couplings begin to 
grow rapidly. It can be estimated from an approximate 
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solution of the second ordet-scaling equation (C9).^^ The 
result found for T^ by vzll3'Pi53'^'^«(''") is 



T,. = i5K(0K(0)]i/' 



w^(0) 
4v^(0) 



(C16) 



Note that the factor [u^(0)w^(0)]^/^ is absent^"* if one 
estim atea _TK only from the leading-order scaling equa- 
tion ^^,p-'^577,eq.{4.ii)^ gijj(,g ^Yie bare v''{0) < 1, this 

factor causes a substantial suppression of T^ (by about 
two order s of magnitude), if one simply inserts w^(0) 
into Eq. (Clt), leadias-|to pessimistically small values 
of Tk — 0.01 — 0.1 KeZIc3. However, the inclusion of ex- 
cited states remedy this problem, in that excited states 
renormalize v^ to larger values by about two orders of 
magnitude (see below). 



5. The Role of Excited States 

Let us now return to t he m ore general problem where 
the excited states in Eq. (P3^ ) with energies £'e^, are not 
neglected from the beginning. 

The first important consequence of including excited 
states in the model has already been discussed in section 
22: electron-assisted hopping transitions between the two 
wells via excited states allow Kondo physics to occur 
even if the barrier is so large that direct and electron- 
assisted tunneling through the barrier is negligible (i.e. 
Ao ~ 0). This is good news, since the energy splitting 
Af, — El ~ E2 is limited from below by Ao, but simulta- 
neously Ab (being a relevant perturbation) needs to be 
very small if scaling to the 2-channel fixed point is to take 
place. 

Secondly, in the presence of excited states, poor man's 
scaling towards strong-coupling, based on Eq. (|C^), has 
to proceed in several steps: the excited state |^^ only 
contributes as long as the effective bandwidth D' is larger 
than Ec, as is made explicit by the 9{D' — Ec) in Eq. (C2). 
As soon as D' < Ec, the excited state decouples. 

Assuming that the presence of excited states does not 
affect the result found in section |C3| , namely that the 
effective Hamiltonian scales towards a 2-D subspace in 
which the piprtrnns ha.vp pseudo-spin Se — ■k, Zarand 
and ZawadowskieZlcS have analyzed the successive freez- 
ing out of excited states. They concluded that when D' 



becomes smaller than the smallest excited-state energy 
E^ , one ends up with a tunneling center of formally pre- 
cisel y th e same nature as the one discussed in sections C 3 
and C 4 , but with renormalized couplings. 

The renormalized couplings turn out to be still small, 
which means that the perturbative scaling analysis of 
sections C 3 and C 4 still applies; however, v^ and v^ 
are renormalized upward by a factor of up to 50 from 
their bare values (which were three orders of magnitude 
smaller than v^ see Eq. (CI)). This has very i mpor - 
tant consequences for the Kondo temperature Eq. (C16), 
which strongly depends on v^: with realistic choices 
of parameters, the Kondo temperature turns out to be 
about 2 orders of magnitude larger with than without 
excited states in the model, and Kondo temperatures 
in the esQerimentally relevant range of 1 to 3 K were 
obtainedH*"''''^". 

To summarize: the inclusion of excited states in the 
model leads to more favorable estimates of the important 
parameters Aq (can be zero) and T^ (larger); but since 
the excited states eventually decouple for small enough 
effective band- widths, they do not affect the flow toward 
the 2-channel Kondo fixed point in any essential way. 



APPENDIX D: RECENT CRITICISM OF THE 
2-CHANNEL KONDO SCENARIO 



Very recently, the claim that Zawadowksi's non- 
magnetic Kondo problem will renormalize to the non- 
Fermi-liquid fixed point of the 2-channel Kondo model at 
sufficiently low temperatures ha s been called into ques- 
tion in two separate papersEjEj. We ignored the con- 
cerns stated there when introducing our NTKM in sec- 
tion n, because there our attitude was phenomenological 
and our aim merely to write down a phenomenological 
Hamiltonian that accounts for the observed phenomena. 
However, the question as to whether or not the bare non- 
magnetic Kondo model does indeed renormalize toward 
the non-Fermi-liquid fixed point 2CK model is an inter- 
esting theoretical one in its own right, which, in our view, 
has become all the more relevant in the light of the ap- 
parent success of the NTKM in accounting for all aspects 
of RB's experimental results. Therefore, we summarize 
the relevant issues in this appendix. 



■^^ Since Tk is only a statement about the onset of rapid 
growth of coupling constants, the value obtained from scaling 
equations derived by perturbation theory is expected to give 
approximately the correct scale even though the scaling equa- 
tions thmrtselves become invalid when the couplings become 
too largeH. 



# 



■^"^The presence of the prefactor to the exponent in (C16) 
course a well-known feature of second-order scaling, see 

ggf,4e5.(3.47) 
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1. Large A due to Static Impurities 

Wingreen, Altshuler and Meir have recently arguedtSi 
that tunnehng centers with very small splittings (A < 
IK) can not occur at all in a disordered material if the 
TLS-electron coupling has the large values that apply 
to the over-screened 2-channel Kondo fixed point. Their 
argument goes a s fo llows: 

If Htc of Eq. (B2) is truncated to the lowe st tw o states 
and written in the left-right basis [see Eq. (C4)], it has 



the following general form: Htc = ^J2A=x,y,z^A'J , 
where Az is the asymmetry energy and A^; the sponta- 
neous hopping rate (for the bare system, time-reversal 
symmetry enables one to choose Aj, = by choosing real 
eigenf unctions, but under renormalization Ay ^ can be 
generated, see below). Hence A^ can be be interpreted 
as an effective local field at the TLS site (in the language 
of the magnetic Kondo problem, this would be called a 
local "magnetic field"). The energy splitting is of course 
A = {Al + Al + Al)y\ 

Now, WAM pointed out that ordinary elastic scatter- 
ing of electrons off other static defects in the system, 
which causes Friedel osciliations (wavelength l/kp) in the 
electron density (see e.g.E3), will make a random contri- 
bution to A A (not considered in Zawadowski's theory). 
WAM suggested that the magnitude of this effect can be 
characterized by the typical splitting A that arises for a 
typical TLS in the presence of static disorder. 

As estimate for A, WAM used (A^)^/^, the average of 
A over all realizations of disorder. Calculating this us- 
ing simple 2nd-order perturbation theory in the coupling 
between the electrons and static impurities, they found^^ 

Ai^e^v^l^fkJ.. (Dl) 

Here i is the mean free path (a measure of the concen- 
tration of static impurities), and f^ the effective TLS- 
electron coupling strength in Eq. ( BIO ). Moreover, WAM 
argue that because A^ has three components, the prob- 
ability distribution P(A) goes to at A = 0, because the 
probability to simultaneously find all three components 
A^ = is vanishingly small. 

WAM estimated v^ — 0.1 at the 2CK fixed point, W 
using a Kondo temperature of about 4K (as cited inEl) 
in the standard formula Vj^ — Xo'gik^T-^jtp) obtained 
from the leading-order scaling equations. (This value 
for Vn agrees, with the values estimated by Zarand and 
Zawadowskilla for the 2CK fixed point (see caption of 
Fig. H). Using I ~ 3 nm, WAM then found a value of 
A c± lOO i^ (the result is so large because according to 
Eq. ( |D1[ ) a is proportional to e_p). 



Since 100 K is a huge energy scale compared to all other 
scales of interest, WAM argued that the 2-channel Kondo 
physics evoked in paper I to explain the Ralph-Buhrman 
experiment would never occur. Instead, they proposed 
an alternative explanation of the experiment based on 
disorder-enhanced electron interactions. The latter sug- 
gestion, psihich we believe contradicts several experimen- 
tal factsE^^*-^', is critically discussed Appendix A 1 of pa- 
per I. Heie we briefly comment on their estimate of A, 
following C") andEll. 

We believe that WAM are correct in pointing out that 
static disorder interactiond can act to increase the en- 
ergy splitting. A, of the TLS. However, we suggest that 
A ~ 100 K may be a considerable overestimate, for the 
following reason. 

Firstly, it was pointed out by SmolyarenkoEJ that the 
quantity (A^)^/^ calculated by WAM does not corre- 
spond to the typical splitting (i.e. the splitting at which 
the distribution P(A) of splittings peaks), but instead 
to the average splitting JdAP(A)A. This distinction is 
important, since the average of A over all realizations 
of disorder can be dominated by a few rather rare real- 
izations of disorder that give rise to very large splittings 
(e.g. where some static defect is very close to the TLS, so 
that the Friedel oscillations have very large amplitudes at 
the TLS), and hence can be much larger than the typical 
splitting. However, for the experiments at hand, clearly 
the relevant quantity is the typical splitting, since TLS 
that have a very large splitting due to a very nearby static 
defect would simply not exhibit Kondo physics. Smol- 
yarenko was able to estimate the typical splitting, and 
found a result of A ~ 15K, much smaller than WAM's 
value for the average splitting. 

Fucthermore, according to Zarand and Zawadowski 
(ZZ)Lil, WAM's statements are equivalent to assuming 
that Aa is renormalized by Hartree^type corrections to 
the TLS self-energy (see Fig. 2 ofEil for the Feynman 
diagram): ^^(t^) = Jd^Paa'ii^)Va>a[^'n-iuj/D)]. Here 
Paa' (i^) is the spectral function of the conduction elec- 
trons in the presence of impurities and v the renor- 
malized vertex function of Eq. ( p4| ) [with x — \n{uj/D) 
there; the corresponding bare vertex function would be 
zJ'^(O)]. WAM's (and Smolyarenko's) estimate of A is 
obtained if one simply uses the unrenormalized diago- 
nal part of the spectral function. However, this is too 
simplistic, since if the renormalised. spectral function is 
used, Aa is reduced significantly^'c3 (despite the growth 
in the couplings v under renormalization). The spon- 
taneous hopping rates A^^Ay, in particular, decrease 
by as much as three orders of magnitude under renor- 
malization (ZZ estimate their final typical value to be 



^^Cox has reproduced WAM's resultE3 by a simple calcula- 
tion analogous to the one by which one obtains the RKKY 
interaction between two magnetic impurities. 
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A^: ~ Ay S 1 K to 0.1 Ktil). This simply reflects the 
screening of the TLS by conduction electrons: when tun- 
neling between the wells, the tunnling center has to drag 
along its screening clowd, which becomes increasingly 
difficult (due to the orthogonality catastrophy) at lower 
temperatures. Irucpntrast, the asymmetry term Az is not 
reduced as muchE3 (ZZ estimate that after renormaliza- 
tion Az ~ 1 K), because of a much larger value of the bare 
coupling w^(0) [~ 103^^(0) ~ 103w!^(0), see Eq. (CI)]. 

Thus, we believe that the reason for the large estimates 
of A is the neglect of the reduction of A^ under scaling. 
Though ZZ's studies of this reduction were performed 
without considering static disorder, disorder should not 
essentially change matters^^ (since this reduction simply 
reflects the well- understood physics of screening). More- 
over, because A^,, Ay end up being so much smaller than 
Az, WAM's conclusion that P(0) = for the distribution 
of splittings is not persuasive, because the distribtions 
Pa{^a) of the individual Aa are not equivalent, as they 
assumed. 

Note that although ZZ estimated that A^ ~ 1 K, im- 
plying that also A ~ 1 K, this is only a statement about 
the typical splitting of a typical TLS. In a disordered 
system, it seems very likely that some TLSs will exist 
with a splitting A significantly smaller than the typical 
A. In particular, ZZ's estimate that typically A^,, Aj, are 
~ 0.1 K implies that assumption (A2) of paper I, Sec- 
tion V D, namely that the nanoconstriction does contain 
TLS with A < 1 K, does not seem unreasonable, despite 
the fact that A ^ 1 K. 

Finally, note that two-state systems with small energy 
splittings (A < IK) have been directly observed in dis- 
ordered metals in at least two experiments: Graebner et 
alxJ found a linear specific heat in amorphous supercon- 
ductors below Tc in the regime 0.1 < T2 K, which they 
attributed to two-state systems; and Zimmerman et al. 
directlyj-observed individual slow fluctuators in Bismuth 
wiresE3'EJ, with As as small as 0.04 K. Though the de- 
tailed properties of these two-state systems may be dif- 
ferent than those of fast TLS, this does illustrate that 
even in systems where the average splitting is expected 
to be large, the physics can be sometimes dominated by 
those two-state systems that have smaller splittings. 

The relevance of WAM's calculation to the interpreta- 
tion of the experiments discussed in paper I are discussed 
in Section V C 3 of paper I. 



2. Another Relevant Operator 

The theoretical justification for the non-magnetic 
Kondo model proposed by Zawadowski has recently also 
been questioned by Moustakas and Fisher (MF)E3. Reex- 
amining a degenerate two-level system interacting with 
conduction electrons, they argued that the model of 



Eqs. (B6) and (B7) used by Zawadowski is incomplete, 
because it neglects certain subleading terms in the TLS- 
electron interaction that have the same symmetries as 
the leading terms. MF showed that when combined in 
certain ways, these subleading terms generate an extra 
relevant operator, not present in Zawadowski's analysis, 
which in general prevents the system from fiowing to the 
T = fixed point. Therefore, unless a fine-tuning of pa- 
rameters miraculously causes this new relevant operator 
to vanish, it will eventually always become large, and the 
system will never reach the T — fixed point. 

Zawadowski et allEj have recently reinvestigated the 
nature of this new relevant term within the context of 
a somewhat simpler (commutative) model than that of 
MF, and suggested that it arises due to the breaking of 
particle-hole symmetry. They estimate that before renor- 
malization, its prefactor in the bare model is smaller than 
the effective Kondo coupling at the flxed point by a fac- 
tor of 10^^, and still by 10^"^ after poor man's scaling 
renormalization to effective bandwiths of order D' — Aq 
(the spontaneous tunneling rate). Thus, they conclude 
that this effect can probably be neglected in realistic sys- 
tems. Moreover, they questioned the path-integral renor- 
malization group scheme employed by MF, arguing that 
it leads to misleading conclusions about the renormal- 
ization of the physical couplings. At present, this mat- 
ter still seeiris-jcontroversial and seems to deserve further 
investigationlI9. 



APPENDIX E: THE LIMIT OF LARGE 
CHANNEL NUMBER: A' ^ oo 

In this appendix, we perform a check on the CFT cal- 
culation of the backscattering current of section ^, by 
considering the (unphysical) limit of a large number of 
conduction electron channels, i = 1, . . . , fc, with k -^ oo. 
In this limit, the poor man's scaling approach becomes 
exact, even though it is based on perturbation theory. 
The reason for this is that (for the isotropic model) the 
over-screened flxed point occurs when the coupling con- 
stant has the special value v* — j+fe; which — > as 



To check this statement, an extra term, including the ef- 
fects of static disorder, should be added to Zawadowski's 
Hamiltonian, and then a full RG analysis should be performed 
to determine self- consistently how the couplings and the "lo- 
cal field" A A, evolve together under renormalization. 
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k -^ oo. Thus, in this limit one never scales into a 
"strong-coupling" regime, and the perturbative expres- 
sions from which the scaling equations are derived retain 
their validity throughout. Therefore, results from the 
poor man's scaling approach should agree with exact re- 
sults from CFT in the limit k —t oo, which serves as a 
useful check on both methods. 

Thus, in this appendix we consider the fc-channel ver- 
sion of the NTKM of section O, in which the index 
i = l,...,fc, but the interaction still has the form of 
Eq. (D. 



1. Backscattering current 

To begin, we need a perturbative expression for the 
backscattering current A/ (i.e. the negative contribution 
to /, which we defined such that J > if it flows to 
the right) due to the backscattering events of V^5 . We 
use the most naive approach for treating the effects of 
the interaction in a non-equilibrium situation: we sim- 
ply do perturbation theory in Hint according to the rules 
of T = 0, y = quantum mechanics, and insert by 
hand, into all sums over intermediate states, appropri- 
ate non-equilibrium distribution functions that indicate 
with what probabilitiy the corresponding states are oc- 
cupied or empty, X),, Idefa{e) or J^r^ X^f [1 - ^(e)], as 



in the definition (|e|) of T. In Eq. {^ the factor 



in Eq. ( B8 ) . This is the method KondoEJ used when de- 
riving his famous log T for the first time, and when one is 
merely interested in the lowest few orders of perturbation 
theory, it is certainly the simplest approach. 

In this approach. A/ is given quite generally by^^ 



A/ 



KNo 



has been included so that ^J2a' represents an average 
over the initial states of the TLS. Tjf" , , is the general- 



ization of the interaction vertex r°°, of Eq. (B8) to all 
orders of perturbation theory, and depends not only on 
the matrix elements Vff, , but also on the distribution 



functions J{e,rf) (as Eq. ( p38| ) illustrates explicitly). The 
factor K = e^ ^ T^{Q)/h (where r,,(e) is defined below) 
is included in Eq. (El) for dimensional reasons, and the 



dimcnsionlcss constants 6,pj' characterize all those details 
of scattering by the impurity that are energy- independent 
and of a sample-specific, geometrical nature, such as the 
position of the impurity relative to the constriction, etc. 
(compare section III of paper I). 

Now, since V^'^5 is independent of the indices c, cr', 
the same is true for |r°°, ,p. (Although the "internal 
sums" '^^11 involving intermediate states are highly non- 
trivial, because of the presence of cr"-dependent f{s,ri") 
functions, T^,'^^,^, of course does not depend on such a"- 
indices, since they are summed over.) Therefore it follows 
immediately that 



r(e', a' -> e, cr) = T{e, a -> e', a') 



(E4) 



Exploiting eq. (E4) and the (5(e'— e) function in r(£', a' -^ 
s, a), the backscattering current can be brought into the 
suggestive form [compare with Eq. (|53|)] 



A/ 



Kb 



fde'[Me')~h{e')]^Y.^nr,- (E5) 



de de' [br,r,'{l - f,{e))f,.{e')T{e\n' 



e, 



6;,^(l-/<,-(e'))/<x(e)r(£,a^£',a')] , 



where a = L and a' = R, and the backscattering rate 
from cr'-movers to cr-movers is 



Here we have taken &,,,,/ = b for simplicity, and have de- 

^ fined-, the total scattering rate — l-pr for a electron with 
a) (El) ^ S'(s) 

energy e' and discrete quantum numbers rj' by^* 
(¥^ r 2vr 



How is T- ^. 



aia.a t a 



(E3) 



In writing down Eq. (El), the fact that the Fermi func- 
tions do not depend on the indices that appear in the 
sums J2aia a'i'a' ^^^ been exploited to pull them out to 
the front of these sums, which could thus be included 



^, to be calculated exphcitly? For T > T^, 
the leading order logarithmic terms of the perturbation 
series in powers of H^nt can be summed up using the poor 
man's scaling approach, discussed in the next subsection. 
On the other hand, an analysis of the regime around the 
T = 0, V = fixed point requires the use of CFT (see 
paper III). A consistency check between the two methods 
can be performed by taking the limit fc — > (X) in which 
perturbation theory becomes exact, as discussed above. 



^'^ Though the relation of the perturbative expression (El) to 
the non-perturbative ones of section M is not readily apparent, 
note that Eq. ( |E5[ ) has the same form as Eq. (p3). 

■^*The ^ in Eq. (Ra) is needed because the definition (Eg) 
of T~^ contains a sum Y^ , that does not occur in Eq. (|E3D, 

and is = 2, since T""/ is independent of a' . 
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2. Gan's Results for Large Channel Number 



1. The KadanofF-Baym Ansatz for \/ / 



A calculation of r^^(e) for the bulk isotropic fc-channel 
Kondo |modcl, in the limit k — > cxi, has been carried out 
by Gar£3. More specifically, he calculated the imaginary 
part of the electron self-energy, 'E'{u!,D,g), perturba- 
tively^^ to order fc~^ (we cite only the lowest relevant 
terms below). By then using poor man's scaling meth- 
ods, he was able to obtain agreement to order k^^ with 
the exact CFT results for T,' . 

Since Gan considered precisely the interaction Hamil- 
tonian of Eq. (n3) that governs our even channels, we 
can directly use his results. He obtained the following 
expression at T = 0, T^ = 0: 



Y.'{uj,D,Vk) oc 



1-Cl 



t7< 



His perturbative expression for the exponent occuring 
here is a = |(1 — I), which agrees to order k~^ with 
the exact CFT result a = ■^^. 



Since Eq. (E7) was derived using poor man's scaling 
methods, it also holds in thenon- equilibrium case, as 
long SiS uj > V (see section |IIE|). This condition does 
not hold strictly in the integra l ( l ESh . Nevertheless, if we 
use -^ = -2E^(w) in Eq. (psT^ith T ^ 0, V ^ 0, 
the resulting asymptotic expression for the backscatter- 
ing conductance AG{V,0) = 9yAJ(F,0), 



AG{V, 0) (X F" , 



should still be approximately correct up to logarithmic 
corrections that are typical of the poor man's scaling ap- 
proach. Indeed, the corresponding expression that we ob- 
tained in section M from our CFT approach has the same 
asymptotic form [see Eq. (pq)], but with a replaced by 
the exact value for the exponent, namely a = jXfe- [Ac- 
tually, in section |V| we always use fc = 2, and hence a — ^ 
in Eq. ( pq ) .] This agreement is a reas suring confirmation 
that the methods used in sections HI to ^ agree with the 
present perturbative results in the one limit (/c — > oo) 
where perturbation theory can be trusted. 



APPENDIX F: HERSHFIELD'S y-OPERATOR 

APPROACH TO NON-EQUILIBRIUM 

PROBLEMS 

In this appendix we summarize the main ideas of Her- 
shfield's y-operator approach to non-equilibrium prob- 
lems. 



The problem at hand is defined by the free Hamiltonian 
Ho of Eq. (0) , the free density matrix po of Eq. (0) and 
the interaction Hint of Eq. (g|). Ho and po describe free 
electrons that move between two leads or baths (R/L), 
at different chemical potentials (y^i), by passing ballis- 
tically through a nanoconstriction, in which L- and R- 
movers can be scattered into each other by Hint- 

How does one calculate statistical averages for such 
a system, in other words, how does one define the full 
density matrix in the presence of Hint"! The main com- 
plication that has to be confronted is that the num- 
ber of electrons in each bath is not conserved, in that 
(E7) [J^L,R,Hint] ^ (compare footnote 6 on pageg). There- 
fore, any attempt to naively replace PoiV) in Eq. (0) by 
Q-0{H-Yo) .^jjj (apart from lacking first-principles justi- 
fication) quickly run into problems: since [H, Yg] ^ 0, 
many of the standard properties of equilibrium Green's 
functions [e.g. G{t + f3) = ±G{t)], no longer hold. 

Kadanoff and Baym have shown how such a general 
problem can be dealt with, by using the notion of adi- 
abatically switching on the interactions'^^' ^^■^*'-': Ther- 
mal weighting has to be done with the initial density 
matrix po at some early time to — > — oo, at which 
all interactions Hint are switched off, and then Hint 
is adiabatically turned on [Hintit) = i/i„ie"*, with 
a — > 0^] while the system is time-evolved to the time 
t of interest. Concretely, to evaluate the thermal ex- 
pectation of an operator O, one writes the operator in 
the Schrodinger picture, and uses the thermal weight- 
ing factors e~'^[-^°"~2^^(^i-^«)"l appropriate to a trace 
X]n("''^ol l^i^o) of Schrodinger states taken at some 
early time to ^ —oo (where they are eigenstates of Ho 
with eigenvalues Eon)- However, one then takes the ac- 
tual trace between the time- evolved versions of these 

states \n,t) — U{t,to)\n,to), whereC/ — Te •'*'' 
is the Heisenberg time-evolution operator: 



(E8) 



(Oit)). 



Er 



,-/3[B„-ieV(JVi,-JVB)„], 



n,t\0\n,t) 



V (>-PlEa,^-^eV{NL-NIi)r,] 

Ti po{V,to)UHt,to)OU{t,to) 



TTPo{V,to) 



(Fl) 



where in the second equality the trace is taken between 
the states \n,to)- (Since steady-state expectation values 
of a single operator are time- independent, t is here just 
a dummy variable, and is often taken to be 0.) 



■^^Since the coupling constant v ~ 1/fe, and closed electron 
loops get a factor fc, Gan had to include up to 8-th order 



diagrams! 
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Eq. ( piD is the defining prescription for taking non- 
equilibrium expectation values in the presence of interac- 
tions. For V^ = 0, it reduces to the standard equilibrium 
prescription,'^*' 



_ Trp(0,OO 
^^^^-°- Trp(0,io) 



where p(0, to) = 



inli. 



~PH 



(F2) 



as shown, e.g., by Hershfield inli3. Eq. (F2) is of 
course the starting point for familiar equilibrium sta- 
tistical mechanics. One of its most useful features is 
that the thermal weighting factor e~^^ and the dynami- 



-iH(t-t„) 



commute; 



cal time-evolution factor U(t,to) 
Green's functions therefore have the periodicity property 
G{t + P) — zLG{t), which makes it convenient to formu- 
late perturbation expansions in i?i„f along the negative 
imaginary axis, t = —ir G [0, —if3]. 



2. Hershfield's Formulation of the case V 7^ 



If V" 7^ so that Eq. (^) and not Eq. (|j) is the 
starting point, there are no obvious periodicity proper- 
ties along the imaginary time axis, and the conventional 
approach, due to Keldysh, is to formulatfi perturbation 
expansions in Hint along the real axisOo. The various 
non-equilibrium diagrammatic techniques that have been 
devised are simply ways of doing the real-time integrals 
/j dt' that result from the expansion of U{t,to)- How- 
ever, for our purposes such expansions are inconvenient: 
firstly, perturbation expansions have limited use in the 
Kondo problem, and secondly, we would in the end like 
to apply Affleck and Ludwig's non-perturbative GET re- 
sults. 



Hershfield has recently shown that Eq. (Fl) can be 
rewritten in a way that exactly meets our needs. The first 
step is trivial: using the cyclica l pr operty of the trace to 
move U{t, to) to the front, Eq. (Fl) can be written as 



(0(i)>v 



Trp{V,t)0 

TTp{V,t) 



where 



p{V,t) _ U{t,to)po{V,to)Wit,to) 



TTp{V,t) 



TTpoiV,to) 



(F3) 



(F4) 



The formal definition (F4) makes it clear that p{V,t) is 
the density operator that Po{V,to) develops into as the 
interaction is switched on and the system time-evolves 
from to to i, with appropriately changing normalization. 
Thus, all complications introduced through the adiabatic 
switch-on procedure are lumped into the time-evolved 
density operator p(V,t). 

Next, Hershfield transfers these complications to a new 
operator, Y, which he defines by writing p{V,t) in the 
form 



p{V,t) 



^ ^-PlH-Y(V,t)] 



(F5) 



purposefully constructed to resemble the definition of 
Po{V,to) in Eq. (Q). Then he was able to show^^ (and 
herein lies the hard work) that the operator Y thus de- 
fined can be characterized as follows: 



(P) 



Y is the operator into which Yq evolves as the inter- 
actions are turned on [as is suggested by a compar- 
ison of Eqs. (f1) and (0)]. It satisfies the relation 



[Y, H] = ia{Yo - Y) 



where 



0^ 



(F6) 



which implies that Y is a conserved quantity. 

The fact that the Y -operator is a conserved quantity is 
the great advantage of the Y-operator approach. It im- 
plies that the problem is now formally equivalent to an 
equilibrium one (for which one has fiN {N— total elec- 
tron number) instead of Y, and [H, N] — 0). 

If Y is known, one can therefore apply the usual 
methods of equilibrium statistical mechanics, ^^ using the 



^"The second argument to in p(0, to) is superfluous; it is re- 
tained here only for the sake of notational consistency with 
the V j^ case. 

^'^Hershfield's proof is perturbative: using Eq. (Fq) he 
showed explicitly that Eq. (F5), expanded in powers of 
Hint, reproduces the Keldysh perturbation expansion ob- 
tained from the KadanofF-Baym Ansatz ([F1|). 

■^^From Eq. (hj) it is clear that Y can actually be shifted 



away m p = e 



r(H-Y) 



by defining new energies e' = e 



Pr, 



associated with Cg,,, i.e. measuring the energy of an excita- 
tion relative to the Fermi surface of the bath from which it 
originates. The weighting factor then completely resembles 



its equilibrium form, but because Csriir) = c^r^e 



r{e'+/J>)) 



tra factors of g*^^'''' appear on some operators that are not 
diagonal in a, such as Hint. We shall not follow this approach 
here. 
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-P(H-Y) 



-,Hi 



Oe 



and Heisenberg time- 
^'^, to calculate physical 



density matrix p 
development 0(r) 
quantities. 

In general, finding Y explicitly is just as difficult as 
solving the full non-equilibrium problem. However, for 
scattering problems, it is sometimes possible to write 
dowUpX explicitly in terms of the problem's scattering 
statesED. Explicit expressions for H a nd Y in a typical 
scattering problem are given in section III A . 



|T|2 + |7e|2 = 1]. Comparing Eq. (|g|) with Eq. (pi) and 



APPENDIX G: EXAMPLE: 2-SPECIES 
POTENTIAL SCATTERING 

In this appendix we illustrate the formalism developed 
in section |IVD by applying it to a very simple scatter- 
ing problem, namely the scattering of only two species of 
(spinless) electrons off a static scattering potential. We 
take T] equal to the species index, rj = a = (L, R) = 
(-1-,—) (i.e. rj contains no extra channel indices «, and 
L/R denotes physical L/R movers). As Hamiltonian we 
take [compare Eqs. (|2|) and (|^]: 

H = Ho + Hint 

= E/ifV^-(*^)('^--'*^- + 2'^'^(^)"^'-')V'a'(«2;') . (Gl) 

Here Vaa' is simply a hermitian 2x2 matrix represent- 
ing potential scattering of the two species into each other 
(i.e. the impurity is not a dynamical object with inter- 
nal degrees of freedom). Since ¥„„' is Hermitian, we can 
make a unitary transformation of the form 

^, = M,,:^„, , (G2) 

with M chosen such that it diagonalizes 

= V^l(O) {vo\5,a- + «3i^L') V^.'(O) . (G3) 

Since the scattering term is now diagonal, its only ef- 
fect on the '0g.-fields can be to cause a phase shift of the 
outgoing fields relative to the incident ones: 






for a; < , (G4) 



where P„„, = 8„ ..-'^i'^o+a^,^) 



and the phase shifts (^q and 03 are functions of Vq and v^- 
Rotated back into the i/'cr-basis, this phase shift of course 
becomes an actual [S'J7(2)] rotation of the two species 
into each other: 

^Ra{ix) = Ucrcr'i>La'{ix) , (G5) 

where t/o-o-' is a unitary matrix of the form: 

U,,, ^ (M-^PM)^^, EE ( _^.^/^. 5 ) , (G6) 



Eq. (^, we see that Uaa'i^') — U^a', i-e. in this sim- 
ple case U is e'-independent. Physically, this rotation of 
physical L- and i?-movers into each other simply refiects 
the fact that Hint causes backscattering: an incoming 
L-mover has amplitude T to undergo forward scattering 
and emerge as a L-mover, and TZ to be backscattered 
into a i?-mover. This illustrates how our formalism is 
able to deal with backscattering despite the fact that we 
expressed both a ~ L and cr = _R as mathematical L- 
movers in Eq. (|2^), for which both the transmitted (T) 
and refiected (7^) parts of ^^ live at x < 0. 

To calculate the current, insert Eq. (G6) into Eq. (|4(|). 
One readily finds 



As expected, the conductance G = dvl — xl"^!^ ^^ ^^" 
duced from its customary value for a single channel in 

2 

the absence of scattering, namely ^, by the transmis- 
sion coefficient squared, |Tp. 

Eq. ( |G7| ) can also be used to illustrate that the con- 
ductance assumes a V/T scaling form if the transmission 
coefficient T is energy dep endent. Assume that for some 
reason the T in Eq. (^^) depends on the energy dis- 
tance from the Fermi surface, and can be expanded as 
|Tp = Ao + [e/ep)Ai + {ejepfAi + ... Then the con- 
ductance G = dyl is readily found to be 



G{V,T) 



Ar, 



A.^ — 



3 



eV 



(G8) 



This has the scaling form G{V, T) = G(0, 0) + BT^T{v), 



(1 



47r2' 



is a universal function, and 



where V{v) 
V = eV/T. 

In the 2CK case, a scaling form for the conductance 
arises in a similar fashion, namely from an energy- 
dependence in the transmission coefficient. The non- 
trivial difference is that there we have |T(e)p = Aq + 
AiTi/^f (e/T) , see section 0. 



APPENDIX H: THE NCA APPROACH 



In section 



VI we compared our results to recent 
numerical calculations by Hettler, Kroha and Hersh- 
field (HKH)ll3, who used the non-crossing-approximation 
(NCA) approach to the Kondo problem. Therefore, a few 
words about their work are in order here. 



32 



1. Anderson model used for NCA 

HKH represent the system by the following infinite- C/, 
SU{2), 2-channel Anderson Hamiltonian in a slave boson 
representation: 

p.a.a,i Oi 

+ Y.^- (-^^ ^^^p-"^ + H-^-) • (HI) 

P,CT,CK,2 



The first term describes conduction electrons in two 
leads, (J = {L,R) = (+,—), separated by a barrier and 
at chemical potentials fig- = fi + a^eV. The electrons are 
labeled by a momentum p, the lead index a, a pseudospin 
index a = (1,2), and their Pauli spin i = (t,i), which 
plays the role of channel or flavor index. [Note that for 
V — there is no need to distinguish between L and 
R electrons, and the cr-index can be absorbed into the 
p-index.] The barrier is assumed to contain an impurity 
level Sd far below the Fermi surface, hybridizing (with 
matrix elements Va, with Vl = Vr for our purposes) 
with the conduction electrons, which can get from one 
lead to the other only by hopping via the impurity level. 
/ and b are slave fermion and slave boson operators, and 
the physical electron operator on the impurity is repre- 
sented by d^j = fl,bi, supplemented by the constraint 

Although this picture of two disconnected leads com- 
municating only via hopping through an impurity level 
docs not directly describe the physical situation of bal- 
listic transport through a hole accompanied by scatter- 
ing off two-level systems, the Hamiltonian (HI) can be 
mapped by a Schrieffer- Wolff transformation onto the 
more physical one [Eq. (^)] introduced in section II A , It 
is therefore in the same universality class and describes 
the same low-energy physics, provided that one identi- 
fies the impurity-induced "tunneling current" /(„„ in the 
HKH model with the impurity-induced backscattering 
current A/ in the actual nanoconstriction. 

HKH calculate the tunneling current. 



Itun{V,T) = I dujAaiuj)[foiuj - eV/2) - fo{uj + 6^/2)] , 

(H2) 

where /o(w) = \/{e^^ + 1)7 by calculating the impurity 
spectral function Ad{u}) using the NCA approximatiorO, 
generalized to V^ 7^ using Kcldysh techniques. 

2. NCA and the limit N,k ^ 00 

The NCA technique was historically developed as a 
self-consistent summation of an infinite set of selected 
diagrams, yielding a set of coupled integral equations 
(the so-called NCA equations) for the self-energies of 



the slave fermions and slave bosons, which have to be 
solved numerically. Hence the NCA has conventionally 
been viewed as an uncontrolled approximation, "mjth no 
small parameter. However, Cox and RuckensteinE£l have 
recently clarified the nature of the approximation by rec- 
ognizing that the NCA equations can be obtained by 
taking a certain large- A^, large fc-limit for a generalized 
SU{N)s X SU{k)f Anderson model. This model has 
preci sely the same form as the equilibrium version of 
Eq. (HI) (with V = and no cr-index), but the spin 
and flavor indices take the values a — 1,...,N and 
i = 1, . . . /c, and the constraint becomes 



N 

E 

a=l 



fiu+J2^lb^ = l 



(H3) 



i=l 



(In the limit £/ < 0, Vcr/|e/| <C 1, this becomes equiva- 
lent to the SuIn)s X SU{k)f Coqbhn-Schrieffer model.) 
They calculated the impurity contribution to the par- 
tition function, Zimp, as a functional integral using an 
approximation reminiscent of a saddle point approxima- 
tion, and showed that the limit A^ — > cxd, fc — > cx) at fixed 
7 = k/N yields precisely the NCA equations for this 
generalized model. [It was pointed out in Ref. |3^), how- 
ever, that the model does not have a true large- A^ saddle 
point, so that 1/A^ corrections cannot be computed sys- 
tematically.] Moreover, the scaling dimensions which the 
"saddle point" yields for the leading spin and flavor op- 
erators, As — 1/(7 + 1) and A/ = 7/(7 -I- 1), turn out 
to agree precisely, for arbitrary N and k (with k > 2^ 
with the exact results derived by Affleck and LudwigEJ 
for the U{l)c X SU{N)s x SU{k)f Kondo model, with 
the impurity spin transforming in an arbitrary SU{N) 
representation. This is useful, since As, for example, de- 
termines the leading critical exponent of the conduction 
electron self-energy [and hence of the conductance in the 
RB experiments], as can be seen directly in the CFT 
approach by doing perturbation theory in the leading ir- 
relevant operator. Moreover, for N = k, these exponents 
are actually independent of A^. 

However, it should be noted that many other proper- 
ties are not the same for A^ = fc = 2 as for A^, A; = 00 at 
fixed k/N = 7. These include for example the temper- 
ature dependence of the entropy and the detailed form 
of the frequency dependence of the T-matrix (which de- 
pends on the amplitudes of subleading terms in an expan- 
sion in the leading irrelevant operator), etc. In particular, 
the Anderson model and the NCA approximation break 
particle hole symmetry, whereas the original 2CK model 
and its CFT solution do not. 

These caveats should be kept in mind when when com- 
paring the NCA results of HKH and our CFT results for 
the universal scaling function F of Eq. (|6^) ; they imply 
that it would be unreasonable to expect "perfect agree- 
ment" . Nevertheless, the NCA does have the distinct ad- 
vantage that when combined with the Keldysh technique, 
it deals with the non-equilibrium aspects of the problem 
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in a more direct way than our CFT approach, and is 
able to go beyound the weakly non-equilibrium regime 
{V ^ Tk) (i.e. it incorporat es so me of the corrections to 
scaling mentioned in section |VC| ). For our purposes, the 
NCA is thus to be regarded as an approximation which 
happens to give correctly the critical exponents near the 
non-Fcrmi-liquid fixed point for some (but not all) quan- 
tities, and which thus interpolates between the T <^ T^ 
and T ^ T,f regimes (in the latter it always works well, 
due to its perturbative origin). 



3. Electron Self-Energy 

One would expect that the most direct comparison be- 
tween CFT and the NCA could be obtained by compar- 
ing the retarded self-energy S"(w) for conduction elec- 
trons at F = 0, calculated from the NCA, with that 
from CFT [essentially the function f(x) of Eqs. (^) 
and (59)]. However, the usefulness of this comparison 
is somewhat diminished by the fact that the NCA self- 
energy is not a symmetric function of frequency, which 
is a result of using the asymmetric Anderson model. 
This asymmetry disappears when calculating the con- 
ductance, because Itun{V) = Itun{—V) in Eq. (H2) even 
if Adieu) 7^ Ad{uj), meaning that the zero-bias conduc- 
tance is the more meaningful quantity to compare (see 
next section). Nevertheless, for w < 0, the CFT and 
NCA results agree very well, see Fig. 4 of Ref. |l^. 



+Vk Y. {flfp~^-pQ/N)cl 



pPiCpai 



(H4) 



4. Impurity Spectral Function Ad{ui) 

The NCA result for the impurity spectral function 
Ad{u}) is shown in Figure 0, which is very instructive, 
in that it illustrates what when eV ^ T^ (a regime not 
accessible to CFTl, the Kondo resonance splits into two 
(as also found iio for a related model). (T^ is defined 
as the width at half maximum of the T^ = impurity 
spectral function at the lowest calculated T.) However, 
note that even for V ~Tk, this splitting has not yet set 
in, illustrating that non-equilibrium effects are not im- 
portant for eV < Tfc- This is the m ain ju stification for 
the the approach followed in section IIIB of neglecting 
all V/Tk corrections to the scattering amphtudes. 



5. Possible Improvements on the NCA 

Very recentljv, Parcollet, Georges, Kotliar and Sen- 
gupta (PGKS)E3 have devised a novel approximation 
scheme which avoids the problems of the NCA mentioned 
and is based on yet another generalization 



in Section H 2 



of the multi-channel Kondo model: 



HpGKS - Y £pCpaiCpai 



p,a..f3,i 



Here /„ , a — 1, . . . ,N are a set of Abrikosov fermions 
satisfying the constraint 



N 



/ ^ JaJoi y J 



(H5) 



which describe an impurity spin transforming as an anti- 
symmetric representation of SU{N) of rank Q; the con- 
duction electrons with spin index a = I, . . . , N and flavor 
index i = 1, . . . , A; transform under the fundamental rep- 
resentation of the SU{N)s and 5C/(fc)/ groups. PGKS 
showed that in the limit N — > oo, k —^ oo, Q — > oo, 
with k/N — 7 and qq — Q/N held fixed, the model has 
a true saddle point, and that its saddle point equations 
are again identical with the NCA equations (bosonic op- 
erators bi corresponding to those of the NCA arise via a 
Hubbard-Stratonovitch decoupling of the Kondo interac- 
tion). Nevertheless the approach of PGKS, though rem- 
iniscent of the NC A, is not fully equivalent to it, since 
the constraint (H5) evidently differs from the constraint 
( |f^ ) of the NCA. [For Q = 1 the constraints would be 
equal if in the NCA the limit e^ — > — oo is taken, since 
then the boson states are energetically so unfavorable 
that they are never occupied, implying that the second 
term in (H3) can be dropped. In this limit, the NCA 
becomes particle-hole symmetric] 

The scheme of PGKS shares with the NCA the prop- 
erty that in the limit iV — > oo, fc — > oo, the scaling expo- 
nents As and Af for the spin and flavor fields are equal 
to the exact values -wtt and Trrr of CFT. In addition. 



N+k 



the PGKS scheme has several advantages over the NCA. 
Firstly, for qq — 1/2, particle-hole symmetry is main- 
tained for all iV and 7, including the saddle point at 
N = 00 and the 2-channel Kondo case of present in- 
terest, N = k — 2, Q = I, for which (H4) is just is the 
standard 2-channel Kondo Hamiltonian in the Abrikosov- 
fermion representation. Secondly, since it is based on a 
true saddle point, 1/A^ corrections can systematically be 
incorporated by calculating fluctuations about the sad- 
dle. Thirdly, if generalized to the non-equilibrium case, 
it could be used to obtain a 1/iV expansion also for 
the non-equilibrium properties of interest in the Ralph- 
Buhrman nanoconstriction experiments, since the non- 
equilibrium conductance can be computed within a La- 
grangian (Keldysh) framework. More precisely, by in- 
cluding fluctuations about the saddle, it should in prin- 
ciple be possible to obtain an expansion in powers of 
1/A^ of the universal scaling function F in Eq. ( |63| ) [this 
expansion should not be confused, though, with that of 
Eq. M)]: 



G{V,T) - G{0,T) 



00 ^ 



\'r,[T/T^,V/T] 



(H6) 



j=o 
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The advantage of such an \/N expansion is that each 
term would go beyond the scahng regime, i.e. cover the 
entire crossover of the Kondo model, from weak all the 
way to strong coupling, or correspondingly from T ^ Tk 
to T ^ Tk- Moreover, in the extreme scaling regime, i.e. 
in the limit T <C Tk, we expect that Fq should coincide 
exactly with the CFT result for F. The reason is that 
the latter is in this limit determined completely by the 
scaling dimension Ag, and that in PGKS's approach this 
dimension is reproduced at the saddle point, at which 
all 1/N corrections vanish and all their results become 
exact. 

It would thus be very interesting to repeat HKH's cal- 
culation using PGKS's approach. 



APPENDIX I: 1/-DEPENDENT CORRECTIONS 

TO U„^, 

In this appendix, T^-dependent corrections to the scat- 
tering amplitudes C/^^/ are discussed. 

A key assumption made throughout this paper was 
that the scattering amplitudes that describe scattering 
in the non- Fermi- liquid regime are V- independent, for 

However, a simple poor 



reasons given in section [II B 



man's scaling argument shows that this assumption can 
not be correct in general: If T^ > T, then the RG flow 
will eventually be cut off at an energy scale of order V. 
In the poor man's RG approach, this is implemented by 
replacing the renormalized bandwidth D' by V in the ef- 
fective interaction vertex. This means that the effective 



renormalized Hamiltonian now is explicitly ^-dependent, 
implying that the same will be true for its scattering am- 
plitudes. 

Intuitively, the F-dependence arises because when V ^ 
0, the difference in Fermi energies of the L- and R 
leads-fjauses the Kondo peak in the density of states to 
splitBEj into two separate peaks (at-energies /i ± ^eV, 
(see Fig. ^ of Appendix^, taken fromtS). However, in the 
non-Fermi-liquid regime, this ^-dependence can never- 
theless be neglected, because when V/Tk ^ 1, the split- 
ting of the Kondo peak by eV is negligible compared to 
its width, which is oc T^ (said differently, then V ^ 
cuts off the RG flow at a point sufhciently close to the 
non-Fcrmi-liquid fixed point that the latter still governs 
the physics). 

To investigate the onset of Kondo peak-splitting ef- 
fects with increasing V but still in the non-Fermi-liquid 
regime, we use the same kind of arguments as the ones 
used by AL to find the leading T/T^ term in G,,,,' (see 
paper HI): V ^ breaks a symmetry of the system 
(namely a = L -^^ R), and the breaking of a symme- 
try allows boundary operators to appear in the action 
describing the neighborhood of the fixed point that had 
been previously fpsbidden (for extensive applications of 
this principle, seeEl'^^'=**°"^^^'^). 

To find the form of the leading V ^ boundary op- 
erator, we argue as follows: V enters the formalism only 
via Yo [see Eq. (||)], which takes the following form when 
written in terms of the fields tpri of Eq. (gO|) or f/;^ of 

Eq. My. 



/oo 



(11) 
(12) 



[For the second line we used (TVcr^iV^^)^^' = (J%^,, 
with N given by Eq. ([lO|).] Eq. (||) shows that Yq mixes 
even and odd channels. Since the CFT solution was for- 
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mulated only in the even sector, the present niodel^'^ can 
strictly speaking not be solved exactly by CFT for V y^ 0, 
unless one neglects all effects of Yq on the fixed point 
physics. 



V 



The form (|^) for Yq leads us to conjecture that the 
leading boundary operator appearing in the action when 
1^ 7^ will have the form^'* 



^Sv ^^WTfrYl 



T^ 



dr 



^eai{T)^oodi'r)) + V'oa»(T)V'ea»(T) 



J K 



drJeoir) 



(13) 



Since the "even-odd" current Jeo defined on the right- 
hand-side has scaling dimension 1 (i.e. ai = l,ao = 0), 
dSv has scaling dimension zero (see paper III), and is 
therefore a marginal perturbation. This means that even 
for V/Tk <C 1, if T/V is made sufficiently small, the sys- 
tem will eventually flow away from the non-Fermi-liquid 
fixed point, at a cross-over temperature T*, say. How- 
ever, since this perturbation is marginal, it only grows 
logarithmically slowly as T is decreased, so that T* will 
be very small. Therefore, the non-Fermi-liquid regime, 
in which one has both V,T -^ Tj^ and T > T^, can be 
rather large. The lack of deviations from scaling in the 



data for the low-T regime (see Section VI of paper I) in- 
dicate that T* is smaller than the lowest temperatures 
obtained in the experiment. 

How does SSv affect the scattering amplitudes? First 
note that the T^-dependence of SSv enters in a very 
simple way, namely as a "parameter" that governs the 
strength of the perturbation. Therefore, the methods of 
section IVD, which extract C/,p,(e) from an equilibrium 



CFT Green's function, are still applicable. 

Standard CFT arguments show that the effect of ^^v 
on Gfjfj' is to simply cause a rotation'^^ of the a = e/o 
indices of the outgoing ry-fields relative to the incident 



However, related models exist which can be treated ex- 
actly by CFT even if V t^, for example the model used by 
Schiller and Hershfield inll3. There, the pseudospin index is 
also the L-R index (i.e. the interaction matrix elements are 



5"), which means that 



Yo 



h^vT. 



i^iii-i'X)a'^^>ipc'iiix) = eV 



Js {ix) 



where J^ is the spin current. Now, in this case it is easy to 
find the exact Y-operator in the presence of the Kondo in- 
teraction. Y must both commute with H and reduce to Yo 
when the interaction is switched off. This is evidently sat- 
isfied \yy Y — eV J_ ^Jsl, where ^7/^, is the z-component 

of the new spin current Jsiix) = Jsiix) + 2n5{x)S (see pa- 
per III). Thus, in the combination H — Y that occurs in the 
density matrix p, eV simply plays the role of a bulk magnetic 
field in the z-direction,rwhich can be gauged away exactly by 
a gauge transformationEj''^''' ^•^■■^'''> , Hence in this model, non- 
equilibrium properties can be calculated exactly using CFT. 

^"^It is easy to check that the operator Jeo is indeed allowed 
at the boundary: it must be the product ^e'&o of boundary 
operators in the even and odd sectors, with quantum numbers 
(charge,spin,flavor)= (Qejeje) = {-Qojojo) = (±1, i, i). 
"l>o, which lives on a free boundary (since the odd sector is 
free), is simply the free fermion field ipoai in the odd sector. 
$e must live on a Kondo boundary, which indeed does allow 
a boundary operators with the quantum quantum numbers 
(±1, |, i), as may b«, checked by AL's double fusion proce- 
dure (see table Ic ofta). 

^^See, for example,ll3. At T = 0, one can prove that SSv 
generates such a rotation by closing the J_ dr integral along 
an infinite semi-circle in the lower half-plane (this-^s allowed, 
because Jeo{z) ~ z~^ -^ along such a contouiO'^''"'^^'); 
having closed the contour, SSv has precicely the form required 
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7]'-fields by 



^f;f;' {V) = 6aa' fe 



cos6'v — isin( 
-z sin 01/ cos 01 



(14) 



Oy = arctan I ;^ I ■ 

Here ^vi-Jf simply a convenient way to parametrize the 
rotationEj, and c is a constant. Thus, the effect of 6Sv 
can be incorporated by replacing the scattering ampli- 
tude C/f^f^/ (e) of Eq. (H) hyRfj^"{V)Uf,"f^'{e). Evidently, 
the final scattering amplitude Uf,n{e) of Eq. ( [45| ) will now 
be y-dependent. 

It turns out that for the simple form (||) used for the 
backscattering matrix Vo-o-' , this extra F-dependence "ac- 
cidentally" cancels out^® in Eq. (|l]) for P^(e), so that 
Eq. (|5l| ) remains valid as written. However, for more 
general forms of V^a' , it survives. To lowest order in 
V/Tji, there will then be a contribution to the conduc- 
tance of the form {V/Tj,)T^/'^Ti{V/T) = T^/'^T2{V/T). 
However, this is evidently only a subleading correction to 
the scaling function of Eq. ( |5^ ) . It is of the same order 
as corrections arising from subleading irrelevant opera- 
tors of the equilibrium theory, that we have argued in 
section V C would not be worth while calculating since 
there are too many independent ones. 

To summarize the results of this appendix: when 
y 7^ 0, corrections to the scattering amplitudes Ujjrj' (e') 
of order V/Tk can arise; however, they only give rise to 
subleading corrections to the scaling function T[v). 
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FIG. 1. A symmetrical square double well potential 
(heavy line), such as that used by Zarand and Zawadowski 
for their model calculations, and the wave-functions for the 
states |r), \l) and the first excited state l^s). 
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FIG. 2. (taken boidB'^'^- ^): The T = electron 

distribution function /i (r) shown (a) at the hole and 
(b) at two points near the hole. The picture is a posi- 
tion-momentum space hybrid, showing the momentum-space 
distribution function /i with its origin drawn at the posi- 
tion r to which it corresponds. A finite temperature simply 
smears out the edges of the two (R/L) Fermi seas. 




2 3 4 
(eV/kBT)V2 



r(r) 



Fermi sea bottom 



-4a -2a 



2a 4a 



FIG. 4. Scaling plots of tltg conductance for (a) the NCA 
calculations of Hettler et a/.li3 and (b) experiment (sample 
#1). With Bs determined from the zero-bias conductance, 
G(0, T) = G(0, Q)+B^T'^^'^ [Eq. (60)], there are no adjustable 
parameters. The temperatures in the NCA- and experimental 
plots are in units of Tk and Kelvin, respectively. 



FIG. 3. (taken from 



5|,F»9. 6-j. rpjjg electrostatic potential 
energy e(^'"'(r), which defines the bottom of the conduction 
band, near a point contact with radius a, shown along the 
2-axis for the case eV > 0. Within a few radii a from the 
hole, e^'"' (r) changes smoothly from —eV/2 on the left to 
+eV/2 on the right. 
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FIG. 6. The conductance scaling function r(i;). Curves 
1,2,3 are the experimental curves of Fig. 11(b) of paper I. 
Curve 4 is the CFT prediction from Eq. (p6|). Curves 5 and 
6 are the NCA results of HKH, with T/Tk = 0.003 and 0.08, 
respectively. All curves have been rescaled in accordance with 
Eq. (|§). 
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FIG. 5. Comparison between NCA theory and experi- 
ment for three individual conductance curves from sample 
# 1. By using Tk as a single fitting parameter and choosing 
Tk = 9>K for sample 1, this kind of agreement is achieved 
simultaneously for a significant mwnber of individual curves 
[Hettler, private communication] ,E3. The NCA curves shown 
here correspond to T = G.2,Tk = 2.4K, 0.2Tk = 1.6K and 
O.lST/f = 1.2K (NCA curves for the actual experimental tem- 
peratures of T = 2.257K, 1.745K and I.IK were not calcu- 
lated.) 
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FIG. 7. (a) The second-order vertex corrections that con- 
tribute to Eq. (B8) and generate the leading order scaling 
equation (C3). (b) The impurity self-energy correction and 
the third-order next-to-leading-logarithmic vertex correction 
that generate the subleading terms in the second-order scaling 
equation (CSf). (Note that subleading diagrams that are gen- 
erated by the leading-order scaling relation derived from the 
diagrams in (a) have to be omitted.) Dashed and solid lines 
denote impurity and electron Green's functions, respectively. 
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FIG. 8. 
a^ = Tr[(v'*) 

case Nf = 3. All th ree norms-itend to the same value, in 
accord with eq. (C12). ConsultcS, from which this figure was 



Scaling trajectories of the matrix norms 
""^^ {A = x,y,z), calculated numerically for the 



taken, for details regarding the initial parameters used. 



FIG. 9. The Kondo resonance in the impurity spectpal 
function Ad{uj), calculated T/Tk = 0.001 using the NCaH. 
For our purposes the most important feature of this figure is 
that the Kondo peak does not start to split for eV < Tk- 
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